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ABSTRACT 


The  refraction  of  wave  packets  by  ocean,  surface  currents  is 
investigated  using  a  refraction  method  based  on  Snell's  law  with 
geometric  group  velocity  G'  ®  U'cos<|>  (Breeding,  1978)  where 
U'  »  dcj'/dk,  <J>  is  the  difference  in  the  direction  of  the  wavelets 
and  packets,  and  the  prime  denotes  motion  relative  to  the  current. 

This  refraction  method  requires  curvature  equations  for  the  directions 
of  the  wavelets  and  wave  packets,  and  an  advection  equation  for  the 
direction  of  the  ray.  Three  equations  are  necessary,  since  after 
refraction  the  wavelet  and  packet  directions  differ,  and  a  ray  is  not 
generally  orthogonal  to  the  wave  fronts  in  a  current.  These  equations 
are  derived  and  applied  to  idealized  current  models.  The  current  models 
are  patterned  after  the  Gulf  Stream  and  the  Circumpolar  Current.  Both 
parallel  and  curved  current  contours  are  considered.  From  this  study, 
it  was  found  that  wave  packets  can  be  totally  reflected  or  deflected 
from  their  original  path  by  a  current  with  the  short  period  waves 
being  affected  to  a  greater  extent  than  the  long  period  waves.  The 
curvature  of  a  ray  was  found  to  have  the  same  sign  as  the  vorticity 
of  the  current  and  to  Increase  as  the  vorticity  increased.  Rays 
were  found  to  bend  away  from  the  current  normal  in  a  following  current 


and  Coward  the  current  normal  in  an  opposing  current  becoming  per¬ 
pendicular  to  the  current  speed  contours  if  sufficient  refraction 
occurred.  Wave  refraction  was  found  to  be  the  least  for  rays  with 

small  initial  angles  and  for  currents  where  u  /G.  (maximum  current 
6  max  i 

speed  to  initial  geometric  group  speed)  is  small.  Further,  it  was 
found  that  the  change  in  the  direction  of  a  ray  was  greater  than  the 
change  in  the  direction  of  the  packet,  which  was  greater  than  the 
change  in  the  direction  of  the  wavelets.  The  ray  angle  increased  in  a 
following  current,  and  decreased  in  an  opposing  current.  Critical 
angles  were  computed  for  parallel  current  speed  contours.  From  an 
analysis  of  these  angles,  it  was  found  that  wave  packets  totally 
reflect  before  the  wavelets,  and  at  the  same  point  as  do  the  rays. 
Also,  critical  angles  were  found  to  be  a  function  of  the  ratio 
u^/G^  (current  speed  at  the  turning  point  to  the  initial  geometric 
group  speed),  and  to  increase  as  this  ratio  decreased.  Further,  it 
was  found  that  rays  with,  initial  angles  less  than  44  degrees  would 
not  be  totally  reflected  by  open  ocean  currents,  since  extremely 
large  current  speeds  would  be  required.  An  analysis  of  rays  in  cur¬ 
rents  with  curved  contours  indicated  that  some  rays  are  focused, 
that  some  rays  exit  and  re-enter  the  current,  and  that  some  rays  may 
be  trapped  within  the  current.  These  are  features  which  cannot  occur 
with  parallel  currents.  Three  additional  refraction  methods  were 
considered  and  compared  to  the  method  presented  here.  They  were  a 
monochromatic  refraction  method  patterned  after  Arthur  (1950),  a  wave 
packet  refraction  method,  using  conventional  group  speed,  developed 
by  Kenyon  (JL971),  and  a  vorticity  equation  used  by  Teague  (1974)  and 


iii 


presented  by  Kenyon  (1971).  It  was  found  that  there  was  very  good 
agreement  in  the  results  of  all  file  methods  for  small  initial  angles 
and  small  current  speeds.  However,  differences  between  the  results 
obtained  by  the  different  methods  increased  greatly  for  mid-range 
angles  (30-75  degrees)  and  for  large  current  speeds. 
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I.  INTRODUCTION 


Ocean  currents  have  been  observed  to  cause  changes  In  the 
direction,  wavelength,  and  amplitude  of  surface  gravity  water  waves 
which  propagate  through  them.  Several  investigations  have  been  made 
in  which  the  interaction  of  currents  and  waves  were  analyzed.  Some 
investigators  used  the  monochromatic  refraction  approach  for  deter¬ 
mining  wave  trajectories.  This  approach  is  based  on  Snell's  law  with 
phase  velocity.  The  phase  speed  is  defined  as  V  «  w/k  ,  where  <u 
is  the  angular  frequency  of  the  wave  and  k  is  its  wave  number. 
Another  approach  which  was  used  determines  the  refraction  of  wave 
groups  and  is  based  on  the  use  of  Snell's  law  with  conventional  group 
velocity  where  the  conventional  group  speed  U  =  dai/dk  is  defined 
as  being  in  the  same  direction  as  the  wavelets  within  the  group.  The 
approach  taken  in  this  study  is  to  determine  the  refraction  of  wave 
packets  using  Snell's  law  with  geometric  group  velocity.  The  geo¬ 
metric  group  speed  is  defined  (Breeding,  1978)  as  G  «  U  cos(0-y)  • 
The  angle  6  is  the  direction  of  the  packet  and  y  is  the  direction 
of  the  wavelets  within  the  packet.  Both  directions  are  defined  with 
respect  to  the  positive  x-axis. 

This  study  will  be  confined  to  idealized,  steady,  non-uniform 
currents  in  deep  water.  Deep  water  is  defined  as  water  of  depth 
greater  than  one-half  the  wavelength  of  the  waves. 


In  Chapter  II,  a  brief  history  of  work  on  wave  refraction  by 
currents  is  presented.  In  Chapter  III,  the  theory  used  in  this  work 
to  determine  the  refraction  of  wave  packets  by  currents  is  developed. 
This  consists  of  curvature  equations  for  wavelets,  packets,  and  rays, 
and  an  advection  equation  for  rays.  Also,  the  physical  properties 
of  rays,  as  well  as  Snell's  law,  and  the  critical  angles  for  the 
wavelets,  packets,  and  rays  are  discussed.  Chapter  IV  contains  des¬ 
criptions  of  the  current  models  used  and  the  trajectories  computed 
for  these  models.  Trajectories  from  parallel  and  non-parallel  current 
models  are  analyzed  and  compared.  Chapter  V  contains  a  comparison 
and  discussion  of  the  results  obtained  in  this  study  and  the  results 
obtained  using  the  refraction  methods  of  previous  investigators. 
Lastly,  Chapter  VI  contains  the  conclusions  and  recommendations. 


II.  BACKGROUND 


Johnson  (1947)  was  one  of  the  first  to  investigate  the 
refraction  of  ocean  waves  by  currents.  He  considered  the  refraction 
of  deep  water  waves  by  a  uniform  current  in  which  the  waves  moved 
with  phase  velocity.  In  Johnson’s  study  wave  fronts  were  constructed 
with  orthogonals  to  the  wave  fronts.  Johnson  showed  that  the 
changes  in  the  wavelength  of  a  wave  refracted  by  a  current  depended 
upon  the  ratio  of  the  current  speed  to  the  wave's  phase  speed  and 
upon  its  angle  of  incidence  to  the  current.  This  relationship  was 
expressed  as 

L/L  ■  (1  -  C  sin  Y  )""2  (1) 

o 

where  L  is  the  wavelength  of  the  refracted  wave,  Lq  is  the  wave¬ 
length  of  the  initial  wave,  C  is  the  ratio  of  the  current  speed 
to  the  phase  speed,  and  y  is  the  angle  the  wave  front  initially 
makes  with  the  boundary  between  the  current  and  still  water.  Johnson 
concluded:  a)  that  the  wavelength,  steepness,  and  direction  of 
travel  of  a  wave  were  significantly  changed  by  its  interaction  with 
a  current,  b)  that  coastal  currents  could  offer  considerable  protec¬ 
tion  against  short  period  waves,  and  c)  that  waves  can  break  within 
the  current  losing  much  of  their  energy  and  thus  creating  an  area  of 
relative  calm  beyond  the  discontinuity. 

Arthur  (1950)  investigated  the  combined  effects  of  current 
and  bottom  topography.  He  considered  the  refraction  of  ocean  waves 
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by  currents  to  be  similar  to  the  determination  of  the  minimal  flight 
path  of  an  airplane.  He  adapted  an  equation  derived  by  Zermelo 
(1931)  for  minimal  flight  paths.  The  air  speed  in  Zermelo*s  equa¬ 
tion,  which  would  be  analogous  to  the  phase  speed  in  wave  refraction, 
was  constant.  Arthur  extended  this  derivation  to  account  for  a 
variation  in  phase  speed.  He  obtained  the  result 


^  3(V'  +  u^)  cosy  3(V*  +  u^)  siny 

3y  3x 


(2) 


where  x  and  y  denote  a  Cartesian  coordinate  system,  y  is  the 
angle  the  wave  direction  makes  with  the  positive  x-axis,  V'  is  the 
phase  speed  of  the  wave  relative  to  the  current,  k  =  2ir/L  is  the 
wave  number,  and  u^  *  (k/k) • u  is  the  component  of  the  current 
velocity  u  in  the  direction  of  V'.  Two  additional  equations  were 
required  to  determine  the  wave  path  Cray)  since,  as  Arthur  pointed 
out,  rays  are  generally  not  orthogonal  to  wave  fronts  in  the  presence 
of  a  current.  These  equations  are 


dx/dt  ■  u  +  V’cosy 

(3) 

dy/dt  ■  Uy  +  V'siny 

where  u  and  u  are  the  x-  and  y-components,  respectively,  of 
x  y 

the  current  velocity.  Arthur  tested  these  equations  on  a  model  of 
a  rip  current.  He  found  that  waves  which  propagated  against  the 
current  slowed  down  and  thus  lagged  behind  waves  of  the  same  period 
propagating  in  the  surrounding  still  water.  He  also  showed  that 
caustics  were  formed  within  the  current.  These  results  were  similar 


to  those  features  observed  for  strong  rip  currents 


Kenyon  (1971)  investigated  the  refraction  of  wave  groups  by 
currents  using  the  geometrical  optics  approximation.  He  determined 
rays  by  integrating  the  ray  equations  (Landau  and  Lifshitz,  1959) 

dk/dt  =  -Vu> 

(4) 

dr/dt  ■ 

for  certain  boundary  conditions.  The  vector  r  is  the  position 
vector  with  respect  to  a  fixed  coordinate  system,  u>  is  the  angu¬ 
lar  frequency  of  a  monochromatic  wave,  V  **  (9/3x)t  +  (3/3y)j  , 

Vk  -  (3/3kx)i  +  (3/3ky)j  ,  and  i  and  j  are  unit  vectors  in  the 
positive  x  and  y  directions,  respectively.  He  also  discussed  the 
physical  implications  of  an  approximate  radius  of  curvature  equa¬ 
tion  developed  by  Landau  and  Lifshitz  (1959)  from  equations  (4)  for 
sound  waves  propagating  in  a  moving  medium.  In  this  equation,  as 
stated  by  Kenyon  (1971),  the  radius  of  curvature  of  a  ray  is 
R  -  U’/C  where  U'  is  the  conventional  group  speed  relative  to  the 
current  and  C  is  the  vertical  component  of  the  current  vorticity. 
This  equation  was  not  used  by  Kenyon  (personal  communication)  to 
determine  rays,  but  was  used  to  discuss  general  properties  of  refrac¬ 
tion  by  currents.  Kenyon,  however,  did  state  that  this  equation  would 
be  a  good  approximate  solution  for  ray  paths  when  the  maximum  current 
speed  was  less  than  5  percent  of  the  initial  group  speed. 

Kenyon  (1971)  used  Gulf  Stream  and  Circumpolar  Current  models 
to  investigate  the  effects  of  refraction  by  a  current.  Kenyon  con¬ 
cluded:  a)  waves  could  be  reflected  or  deflected  from  their  original 
paths  by  a  current,  b)  rays,  for  waves  which  propagate  with  (against) 


a  variable  current,  bend  in  the  direction  of  decreasing  (increasing) 
current  speed,  and  c)  waves  can  be  trapped  about  a  local  minimum 
or  maximum  in  current  speed.  He  also  showed  graphically  what 
Arthur  (1950)  had  implied:  rays  are  not  perpendicular  to  wave  fronts 
when  in  a  current. 

Teague  (1974)  applied  the  radius  of  curvature  expression 
developed  by  Landau  and  Lifshitz  (1959)  to  wave  groups  propagating 
through  a  cyclonic  eddy  and  a  loop  current.  Teague  showed  that 
eddies  and  eddy  type  currents  can  have  the  effect  of  a  giant  lens 
focusing  waves  into  caustics  either  within  or  beyond  the  boundary 
of  the  eddy. 

Hayes  (1980)  applied  Johnson’s  (1947)  equation  for  the 
refraction  of  monochromatic  waves  by  a  current  to  data  obtained  for 
the  Gulf  Stream  from  the  Marineland,  Florida  experiment  (Shemdin 
et  al.,  1975).  From  the  spectral  wave  data  obtained,  waves  with 
dominant  frequencies  of  0.12  Hz.  and  0.078  Hz.  were  observed  near 
the  western  boundary  of  the  Gulf  Stream,  and  a  dominant  frequency 
of  0.12  Hz.  was  observed  for  waves  in  the  shelf  area  west  of  the 
stream.  In  his  monochromatic  refraction  model  the  0.078  Hz.  wave 
was  reflected  and  the  0.12  Hz.  wave  passed  through  the  current  under¬ 
going  only  slight  refraction.  These  results  were  consistent  with 
the  data.  However,  an  analysis  based  on  the  refraction  of  wave 
packets  would  also  agree  with  the  data.  The  0.078  Hz.  wave  had  an 
initial  angle  which  exceeded  the  critical  angle  for  either  the  mono¬ 
chromatic  waves  or  the  wave  packets.  The  0.12  Hz.  wave  approached 
the  current  at  near  normal  incidence.  For  such  angles  of  incidence. 


Chere  is  little  difference  between  the  monochromatic  refraction 
appraoch  and  the  wave  packet  refraction  approach.  Thus,  on  the 
basis  of  Hayes'  (1980)  investigation  it  is  not  possible  to  deter¬ 
mine  the  validity  of  either  the  monochromatic  or  the  wave  packet 
refraction  methods. 

When  there  is  a  variation  of  water  depth,  the  traditional 
method  of  computing  wave  trajectories,  either  for  monochromatic 
waves  or  wave  groups,  is  to  compute  the  paths  using  Snell's  law 
with  phase  speed.  The  wave  groups  are  assumed  to  have  the  same 
trajectories  as  the  monochromatic  waves.  Breeding  (1972,  1978, 

1981)  and  Black  (1979)  have  shown  this  traditional  method  to  be 
incorrect  when  applied  to  the  refraction  of  ocean  waves.  Breeding 
(1972  and  1978)  used  computed  backtracks  from  directional  wave 
data  due  to  Hurricane  Betsy  (1965)  and  Hurricane  Fifi  (1974). 

Black  used  directional  wave  energy  spectra  determined  from  data 
obtained  in  the  Gulf  of  Mexico  due  to  Hurricane  Eloise  (1975).  Both 
Breeding  and  Black  found  that  the  wave  trajectories  determined  by  a 
monochromatic  refraction  law  were  inconsistent  with  the  data. 

Breeding  (1978)  has  proposed  a  refraction  method  for  wave 
packets  (wave  groups)  that  has  yielded  results  consistent  with  the 
data.  He  found  that  after  refraction  the  packets  and  the  wavelets 
within  the  packets  no  longer  move  in  the  same  direction.  This  raises 
the  question  as  to  which  refraction  law  should  be  used  to  determine 
the  ray  paths.  Breeding  (1978  and  1981)  concluded  that  the 
refraction  law  given  by  Snell's  law  with  geometric  group  velocity 
should  be  used  in  computing  ray  trajectories.  The  use  of  this 
refraction  law  has  led  to  backtracks  consistent  with  the  directional 


wave  data  and  the  known  movements  of  Hurricanes  Betsy  (1965)  and 
Fifi  (1974).  The  geometric  group  velocity  G  is  used  instead  of 
the  conventional  group  velocity  U  ,  since  the  wavelet  and  packet 
directions  differ  and  U  is  defined  as  being  in  the  same  direc¬ 
tion  as  the  wavelets.  The  geometric  group  velocity  reduces  to 
the  conventional  group  velocity  when  the  wavelets  and  the  packets 
are  in  the  same  direction.  In  this  study.  Breeding’s  method  of 
wave  packet  refraction  will  be  extended  to  apply  to  the  refraction 
of  wave  packets  by  currents. 


III.  THEORY 


A.  Curvature  equations  for  wavelets  and  packets 

The  refraction  of  wave  packets  by  a  steadily  moving,  variable 
deep  water  current  can  be  computed  by  using  curvature  expressions 
for  the  wavelets  and  packets,  and  an  advection  expression  for  the 
ray.  An  expression  for  the  ray  by  itself  is  not  sufficient  to 
determine  the  ray  path.  The  wavelet  and  packet  directions  along 
the  ray  must  also  be  determined  since  the  trajectory  of  a  wave 
packet  depends  on  the  direction  of  the  wavelets  within  the  packet 
(Breeding,  1978) ,  and  in  the  presence  of  a  current,  a  ray  generally 
is  not  normal  to  a  wavefront  (Arthur,  1950,  and  Kenyon,  1971). 

The  curvature  expression  for  the  wavelets  and  packets  can  be 
derived  from  the  ray  equations  of  geometrical  optics.  The  ray 
equations  for  the  wavelets  and  for  the  packets  are 


dk/dt  *  -Vu) 

(5) 

dm/dt  ■  -7a 

(6) 

where  <*)  is  the  average  angular  frequency  of  the  component  waves 
and  k  is  the  corresponding  wave  number.  The  quantities  m  , 
the  wave  number  of  the  packet,  and  a  ,  the  angular  frequency  of 
the  packet,  are  defined  by  the  incremental  changes  in  k  and  w 
as  (Breeding,  1978) 


It  should  be  noted  that  Ak  involves  the  difference  in  both  magni¬ 
tude  and  direction  of  the  component  wave  numbers.  It  will  be  assumed 
that  the  changes  in  wave  properties  such  as  wavelength,  direction, 
and  amplitude  vary  slowly  over  distances  of  the  order  of  a  wavelength 

The  angular  frequency  relation  for  wavelets  w  *  vk  is  valid 
only  for  waves  in  a  medium  at  rest.  An  angular  frequency  relation¬ 
ship  for  a  moving  medium  was  derived  by  Landau  and  Llfshltz  (1959) 
for  sound  waves,  and  by  Longuet-Higgins  and  Stewart  (1960)  for  ocean 
waves. 

Take  a  rectangular  Cartesian  coordinate  system  denoted  by 
x,  y  fixed  on  the  sea  surface  and  a  similar  system  denoted  by  x',  y* 
moving  with  the  medium  (e.g.,  a  current).  In  the  primed  system, 
where  the  medium  appears  to  be  at  rest,  the  wave  equation  for  the 
packets  has  the  form  (Breeding,  1981) 

ft+€ 

H  *  exp{i(k' •  r'  -  cj't)}  j  A(w)exp{i(m' •  r*  -  a't)}dw  (9) 

Jn-e 

where  n  is  the  vertical  displacement  of  the  sea  surface,  ft  is 
the  average  angular  frequency  of  the  packet,  2e  is  the  frequency 
bandwidth,  and  A  is  the  amplitude  of  the  waves.  In  equation  (9) 
the  quantity  in  front  of  the  integral  represents  the  motion  of  the 
wavelets.  The  quantity  within  the  integral  represents  the  motion 
of  the  packets.  The  position  vector  for  the  primed  system  r*  is 


11 


related  to  the  position  vector  of  the  unprimed  system  r  by 


(10) 


where  u  is  the  velocity  of  the  current.  By  substituting  equation 
(10)  into  equation  (9) ,  the  form  of  rt  in  the  unprimed  system 
becomes 


n  ■ 


f2+e 

(u  •  k'  +  o)f)  t  j  }  |  A(oj)exp{i  £  m'  •  r  - 
ft-e 


(u  •  m*  +  o')  t  J  }du) 
(11) 


Also,  in  the  unprimed  system,  n  must  have  the  form 

ftfe 

n  -  exp{i(k  •  r  -  cut)  j  A(a))exp{i(m  •  r  -  ot)}da)  (12) 

n-e 

Equations  (11)  and  (12)  must  be  equivalent  for  all  time.  In  parti¬ 
cular  when  t  *  0,  It  is  seen  that  |k' |  *  |k|  and  |m'|  **  |m|  . 
Further,  the  angular  frequencies  in  the  fixed  systems  are  seen  to  be 
defined  by 


(D  *  w'  +  u  •  k 

a  -  a’  +  u  •  m 


(13) 

(14) 


Thus,  the  angular  frequency  of  a  wavelet  or  a  packet  in  a  moving 
medium  is  seen  to  be  the  sum  of  the  angular  frequency  of  the  wave 
relative  to  the  medium  and  the  advection  effect  of  the  medium. 

The  curvature  expressions  for  the  wavelets  and  packets 
are  derived  from  equations  (5),  (6),  (13),  and  (14).  The  curvature 


expression  for  the  direction  of  the  wavelets  is  obtained  by  substi¬ 
tuting  equation  (13)  into  equation  (5) .  This  yields 

dk/dt  *  -kW*  -  V(u  •  k)  (15) 

The  second  term  on  the  right  can  be  expanded  using  a  vector 
identity  to  obtain 

dk/dt  =  -kW*  -  k  •  Vu  -  k  x  (V  x  u)  -  u  •  Vk  -  u  x  (V  x  k)  (16) 

The  fourth  and  last  terms  on  the  right-hand  side  of  equation  (16) 
are  zero  since  k  and  r  are  independent  variables.  Hence,  this 
equation  simplifies  to 

dk/dt  -  -kW'  -  k  •  Vu  -  k  x  (V  x  u)  (17) 

An  alternative  expression  for  dk/dt  is 

dk/dt  -  kCd^/dt)  +  e^( dk/dt)  (18) 

A  4- 

where  e^  =  k/k  represents  a  unit  vector  in  the  direction  of  k. 
Equations  (17)  and  (.18)  can  be  combined  and  simplified  to  find 

de^/dt  ■  -W  -  e^  •  Vu  -  x  (V  x  u)  -  e^/k(dk/dt)  (19) 

A  unit  vector  n^  is  taken  tangent  to  the  wavelet  crest,  i.e., 
perpendicular  to  e^.  Thus 

dy/dt  -  •  de^/dt  (20) 

where  y  denotes  the  direction  of  k  with  respect  to  the  positive 
x-axis.  Equation  (19)  can  be  substituted  into  equation  (.20)  to 


obtain 


3(V'  +  u,  ) 


3(V*  +  u,  ) 


2] 


where  V'  is  the  phase  speed  of  the  wavelet  relative  to  the  current 
and  u^  =  u .  This  is  the  same  equation  used  by  Arthur  (1950) 
to  determine  the  paths  of  monochromatic  waves.  This  equation  is 
applicable  to  waves  in  any  depth  of  water.  Arthur  applied  it  to  shallow 
water  waves.  In  this  study  it  will  be  applied  to  deep  water  waves. 

In  deep  water  the  wavelet  speed  can  be  expressed  as 
V*  =  g/ou',  where  g  is  the  acceleration  due  to  gravity.  The 
gradient  of  V*  is  expressed  as 


substituting  equation  (13)  into  (22)  leads  to 

W  =  (3uk/3x)i  +  C3ufc/3y)j  (.23) 

This  value  of  W*  can  be  substituted  into  equation  (21)  to  obtain 


dy  f 3uk  9uk 

=  2  —  siny - -  cosy 


Equation  (2A)  can  be  expressed  in  a  more  convenient  form  by  using 
u^  *  u^cosY  +  u^sinY.  This  leads  to 


3u 

/3u 

3u 

— ^  +  siny  cosy (  — — 

_ 1 

3x 

V3x 

3y 

,  3u 
2  x 
-  cos  Y  - 


Further,  the  arc  length  along  the  trajectory  of  an  orthogonal 


to  a  wavelet  crest  is  defined  by  ds^"  (v'  +  u^Ht 


As  a  result 


equation  (25)  becomes 


dy  2  r  2  3uv  fK  3u 

ds  =  V'  -t-'TT  sin  Y  — ■ 1  +  sinY  cosY - 

Sk  'He  3x  V  3x  3y 


-  3u 
Z  x 
-  cos  y  — — 


This  is  the  curvature  equation  which  will  be  used  to  determine  the 
change  in  the  wavelet  direction  along  a  ray. 

The  curvature  expression  for  the  packets  can  be  obtained 
by  substituting  equation  (14)  into  equation  (6) .  Following  the 
same  procedure  as  for  the  wavelets,  the  result  is 


d0  _  1 

ds  G'  +  u 


3(G*  +  u  )  3(G'  +  u  ) 

sin9  - — - cos0  - — 


where  0  is  the  angle  between  m  and  the  positive  x-axis,  ds 

m 

is  the  incremental  distance  measured  along  the  orthogonals  to  the 
packet  crest,  G'  is  the  geometric  group  velocity  measured  relative 
to  the  current,  and  u^  is  the  component  of  the  current  velocity 
in  the  direction  of  in. 

As  was  equation  (21),  equation  (27)  can  be  rewritten  for  the 


case  of  deep  water  waves.  In  deep  water  the  geometric  group  velocity 
can  be  expressed  as  G'  ■  JsV'cos^.  Making  this  substitution  for 


G*  in  equation  (27)  yields 


ds  G'  +  u 
m  m 


»  3u  (  3u  3u 


sin  9  — L  +  sin9  cos9l - - 

3x  \3x  9y 


2  DU 

-  cos  9  — —  +  —  cos$  .  sin9  (  cosy  — “  +  sinY  — ^ 
3y  1  \  3x  3x 


3u  3u 

x  ,  ...  : 


cosy  -  +  sinY  — ‘ 

3y  3y 


-  COS0 


(28) 


This  is  the  curvature  equation  which  will  be  used  to  determine  the 
change  in  the  packet  direction  along  a  ray. 

B.  Advection  equation  for  rays 

The  determination  of  rays  for  wave  packets  can  be  obtained 
by  taking  Arthur’s  (1950)  approach  to  wave  refraction  by  currents. 

He  recognized  that  refraction  of  surface  waves  by  a  current  was 
due  solely  to  advection.  He  thus  determined  the  wavelet  direction 
from  equation  (2)  and  the  ray  direction  from  the  x-  and  y-components 
of  the  ray  velocity  (equations  3) ,  which  is  the  vector  sum  of  the 
current  velocity  and  the  phase  velocity.  Ray  paths  were  then 
determined  by  simultaneously  solving  these  equations. 

Rays  for  wave  packets  can  be  derived  in  a  similar  manner, 
except  that  both  the  wavelet  and  packet  directions  are  required. 

The  direction  of  the  packets  is  given  by  equation  (28),  and  the 
direction  of  the  wavelets  by  equation  (26).  The  ray  velocity 
expressed  in  its  component  form  is 

dx/dt  ■  u  +  G’cos9 
x 

(29) 

dy/dt  ■  Uy  +  G’sinS 

where  u^  and  u^  are  the  x-  and  y-components  of  the  current 
velocity.  Thus,  the  direction  of  the  ray  p  is  given  by 
.  u  +  G’sin9 

tanp  m  m  —  ■  -  (30) 

“  ux  +  G'cos9 

Ray  paths  for  wave  packets  are  determined  by  simultaneously  solving 
equations  C26) ,  (28),  and  (29). 


C.  Wavelet  and  packet  directions 


Wavelets,  packets,  and  rays  travel  with  different 
velocities  and  after  refraction  in  different  directions.  This 
is  an  important  fact  which  must  be  taken  into  account  when  com¬ 
puting  trajectories.  Even  though  the  directions  of  the  wavelets  and 
packets  differ  from  the  direction  of  the  ray,  they  are  nevertheless 
carried  along  by  the  ray,  and  thus  must  reach  the  same  speed  contour 
as  does  the  ray.  To  compensate  for  these  differences  in  travel 
time  and  direction,  the  curvature  equations  for  wavelets  and  packets 
must  be  multiplied  by  the  ratios  of  the  speeds  and  the  cosines  of 
the  directions. 

Consider  a  current  with  locally,  parallel  speed  contours, 
and  an  x" ,y "-coordinate  system,  such  that  the  x"-axis  is  perpen¬ 
dicular  to  the  contours  (see  Figure  1).  It  is-  thus  seen  that  for 
the  wavelets 

ds^  ■  (P  cosp  dt)/cosy  (31) 

and  for  the  packets 

ds  ■  (P  cosp  dt)/cos0  (32) 

m 

+  ■> 

where  P  *  G'  +  u  is  the  ray  velocity. 

These  equations  relate  the  incremental  distances  traveled 
by  the  wavelets  or  packets  to  that  transversed  by  the  ray.  Sub¬ 
stituting  equation  (31)  into  (26)  and  equation  (32)  into  (28) 
results  in  the  curvature  expressions  needed  for  computations  of 
ray  paths. 


I 

I 
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Figure  1.  Relationship  between  the  incremental  distances  for  the 
ray  (dsr) ,  for  the  wave  packet  (ds,n) ,  and  for  the 
wavelets  (ds^)  between  two  points  along  the  ray.  The 
yH-axis  is  chosen  to  be  parallel  to  the  current  speed 
contours. 


D.  Physical  properties  of  rays 


Some  of  the  physical  properties  of  wave  rays  can  be  deduced 
from  equation  (30).  This  equation  can  be  written  in  a  more  con¬ 
venient  form  by  considering  a  variable,  unidirectional  current  with 
parallel  speed  contours  and  xy-coordinate  system  such  that  the 
y-direction  is  parallel  to  the  current.  Equation  (30)  thus  becomes 

tanp  *  u/(G'cos0)  +  tan0  (33) 

From  this  result  it  is  seen  that  p  >  9  for  a  following 
current  (u  >  0) ,  that  p  <  9  for  an  opposing  current  (u  <  0) ,  and 
that  when  0  *  90  degrees,  the  ray  direction  p  ■  90  degrees.  Thus 
when  the  packets  reflect  the  rays  must  also  reflect.  Also,  it  is 
seen  that  if  u  «  G'  the  ray  and  packet  directions  are  almost 
equal. 

Consideration  of  the  instances  when  p  or  0  are  either 
90  degrees  or  zero  afford  some  interesting  features.  If  0*0  degrees, 
i.e.,  the  packets  are  perpendicular  to  the  current  contours,  the 
direction  of  a  ray  varies  only  as  u/G'  varies.  Thus,  if  the  wavelets 
and  packets  enter  a  current  at  normal  incidence,  they  remain  normal 
to  the  current,  but  the  ray  direction  changes  as  the  current  speed 
changes.  If  the  current  speed  is  constant  the  ray  is  a  straight  line. 
From  the  more  general  equation  (30)  for  p  ■  90  degrees,  either 
u^  ■  G^  ■  0  or  ux  "  -G^,  provided  the  numerator  is  finite  and  does 
not  vanish.  So,  at  a  turning  point  either  the  x-component  of  the 
current  speed  and  group  speed  are  zero  or  they  balance  each  other. 

It  has  been  found  (Figures  14  thru  22)  that  for  currents  with  parallel 


contours,  total  reflection  occurs  for  following  currents  only.  There¬ 
fore,  the  first  case  applies  and  at  a  turning  point  the  directions 
of  the  packet,  current,  and  ray  are  all  parallel  to  the  y-axis.  The 
second  case  will  not  occur  for  parallel  contours,  but  if  the  contours 
are  curved  it  is  possible  for  the  x-components  of  the  speeds  to 
balance  each  other.  In  a  current  with  curved  contours,  it  was  found 
(Figures  26  and  30)  that  a  ray  in  an  opposing  current  can  be  totally 
reflected.  When  the  turning  point  is  reached  under  these  circum¬ 
stances,  the  ray  direction  and  current  direction  exactly  oppose  each 

other.  For  p  ■  0  degrees,  either  u  *  G*  •  0  or  u  =  -G', 

y  y  y  y 

provided  the  denominator  is  finite  and  does  not  vanish.  Only  the 
second  case  is  of  interest.  This  result  indicates  that  rays  may  become 
perpendicular  to  the  current  contours.  It  was  found  (Figures  19  and 
21)  that  for  opposing  currents,  rays  are  turned  toward  the  normal  to 
the  current  and  in  some  instances  become  perpendicular  to  the  current 
contours. 

E.  Snell's  law  for  wavelets  and  packets 

Snell's  law  expressions  for  wavelets  and  packets  are  derived 
by  considering  a  current  with  parallel  speed  contours  such  that  the 
magnitude  of  the  current  velocity  varies  across  the  current.  Each 
current  speed  contour  acts  as  a  discontinuity  to  an  approaching  wave 
causing  the  wave  to  reflect  and  refract.  The  wavelet  or  packet 
must  remain  continuous  across  the  discontinuity.  Therefore,  as  an 
Incident  wave  crest  advances  one  wavelength  the  refracted  wave  crest 
also  advances  one  wavelength.  The  wavelength  changes  as  the  wave 
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crosses  the  discontinuity,  but  the  time,  relative  to  a  fixed  coordi¬ 
nate  system,  required  for  the  incident  wave  and  refracted  wave  to 
advance  one  wavelength  (i.e.,  the  period)  does  not  change.  There 
is,  however,  a  change  in  the  wave  period  measured  relative  to  the 
current. 

The  component  of  the  wave  number  parallel  to  the  boundary 
must  be  constant.  Thus 

'-■*  k  siny  *  constant  (34) 

where  y  is  the  angle  the  wave  crest  makes  with  the  boundary. 

Both  sides  of  equation  (34)  can  be  divided  by  u)  •  Vk.  Therefore 

=  constant  (35) 

where  V  is  the  phase  speed  relative  to  a  fixed  coordinate  system. 
The  phase  speed  V  can  also  be  expressed  as  the  sum  of  the  phase 
speed  relative  to  the  current  and  the  component  of  the  current  in 
the  direction  of  k,  i.e., 

V  •  V'  +  u  siny  (36) 


Then  the  refraction  law  for  the  wavelets  can  be  stated  as 
siny.  siny. 

- 1 - i -  (37) 

+  utsinyt  Vj,  +  u^siny^ 

where  the  subscript  t  denotes  the  transmitted  wave  and  the  subscript 
i  denotes  the  incidence  wave.  If  u^  *  0  this  result  reduces  to 
Johnson’s  (1947)  equation  (3)  which  is 
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siny. 

sinY..  3  - 2  (38) 

(1  -  u  sinY^/V^) 

It  is  clearly  seen  from  this  equation  that  the  refracted  angle  for 
the  wavelets  is  dependent  upon  the  initial  angle  and  the  current 
speed  ratio. 

The  Snell's  law  expression  for  packets  is  obtained  in  a 
similar  manner*  A  crest  which  is  determined  by  the  interference 
maximum  of  the  packet  must  be  continuous  across  a  boundary.  Thus 

m  sinG  ■  constant  (39) 

where  9  is  the  angle  the  packet  crest  makes  with  the  boundary.  In 

a  fixed  coordinate  system,  cr,  the  angular  frequency  of  the  packet, 
is  a  constant  on  both  sides  of  the  boundary.  Thus,  when  both  sides 
of  equation  (39)  are  divided  by  O  the  result  is 

a  constant  (40) 

Further,  G  can  be  expressed  as 

G  ■  G'  +  u  sin0  (41) 


where  G'  is  the  geometric  group  speed  relative  to  the  current  and 
u  sin0  is  the  component  of  the  current  speed  in  the  direction  of  m. 
The  refraction  law  for  the  packets  can  then  be  stated  as 


sin0„  sin04 

C  4. 

-■  1  ■  11  m  —  - 

G£  +  utsin0t  Gj^  +  u^sin9^ 


(42) 


When  there  are  no  currents,  equation  (42)  reduces  to  the  results 
obtained  by  Stonely  (1935)  and  Breeding  (1978).  Stonely  used  the 


result  to  determine  the  direction  of  amplitude  (group)  fronts  along 
monochromatic  trajectories,  whereas  Breeding  used  the  result  to 


define  wave  packet  trajectories.  If  u^  =  0,  equation  (42)  can  be 
rewritten  as 


sin0t 

&t  +  uts*n®t 


sin6 . 

_ i 

G! 


i 


(43) 


Rearrangement  of  this  equation  leads  to 

-  uJ;sinei  G' 
sin0^  sin0t 

Note  that 


(44) 


G'  =  U' (cos0  cosy  +  sin©  siny  ) 
t  t  t  t  t  t 

Therefore 


(45) 


— - —  *»  U'(ctn0  cosy  +  siny  )  (46) 

sin©t  c  t  t  t 

Thus,  substituting  equation  (46)  into  equation  (44)  yields 
G.  -  u  sin©.  -  U'sin0  siny 

ctn0  »  — - - - = - - - - - -  (47) 

U^sin0icosyt 

From  this  equation  it  is  seen  that  the  refracted  angle  for 
the  packets  is  not  solely  determined  from  the  initial  wave  quantities 
period,  angle,  etc.,  but  that  certain  transmitted  wavelet  values  are 
also  needed. 

Equations  (38)  and  (47)  are  useful  in  determining  refracted 


angles  and  critical  angles 


F.  Critical  angles 


A  wavelet  or  a  packet  totally  reflects  when  its  transmitted 
angle  equals  90  degrees.  Thus,  expressions  for  computing  critical 
angles  can  be  obtained  by  solving  the  Snell’s  law  equations  for  the 
corresponding  initial  angles  when  the  transmitted  angle  is  90  degrees 
For  wavelets  when  y  =  90  degrees  equation  (38)  becomes 

siny 

- £ - -  =  i  (48) 

(1  -  usiny  /V  ) 

T  c  i 

where  Yc  is  the  critical  angle  for  the  wavelets  and  u^  is  the 
current  speed  at  the  turning  point.  This  equation  can  be  rearranged 
and  written  as 


The  quadratic  equation  can  be  used  to  solve  this  equation  for  sinYc« 
The  result  is 


siny. 


(2(ux/Vj)  +  l)  -"\/(2(ut/V1)  +  l)2  -  4(ut/V1): 


2(uT/Vi)' 


(50) 


Let  M  ■  u^/G^  "  ^ut^i  ant^  su^st^tute  this  into  equation  (50).  It 
is  found  that 


sinyc  -  ^2(M  +  1)  -  2/2M  +  1  ) /M2  (51) 

Thus,  a  critical  angle  for  wavelets  can  be  computed  given  a  current 
speed  to  Initial  geometric  group  speed  ratio.  Critical  angles  for 
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wavelets  are  independent  of  period  and  shear. 

From  equation  (47)  for  the  packets,  when  0  =  90  degrees, 

ctn0t  =  0  and  it  follows  that 

G.  -  u  sin0  -  U'sin0  siny  =  0  (52) 

i  x  C  T  C  T 

where  =  0^  is  the  critical  angle  for  the  packets.  Then 

G. 

sin0  =  - - -  (53) 

c  u  +  U's iny 
T  T  T 

If  =  0,  (wavelets  and  packets  are  assumed  parallel  upon  entering 
a  current)  then 

sin0  =  U./(u  +  U'siny  )  (54) 

c  l  T  T  T 

This  equation  is  not  in  a  usable  form.  It  is  necessary  to  express 
y^  as  a  function  of  y^  where  =  9c.  Equation  (54)  can  be  written 
in  a  form  independent  of  y^,  if  and  siny^  are  expressed  in 

terms  of  the  initial  wavelet  values.  Equation  (38)  expresses  yt  in 
terms  of  the  initial  wavelet  values.  The  quantity  U|.  can  be  expressed 
in  terms  of  the  initial  values  as  follows 

-  4  /gLT/2tr  (55) 

where  *  L^sinY^/siny^.  Making  this  substitution  for  L^,  and 
equation  (38)  for  siny^,  equation  (55)  becomes 

gLi 

2tt(I  -  u  siny^/V^)^ 

2 

The  quantity  t  ■  gT^/2Tr  and  u/V^  =  ljM  where  denotes  the 
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initial  wave  period.  Thus,  using  these  definitions  for  and  M, 
and  rearranging  terms,  equation  (.56)  becomes 


(1  -  sinyi) 


(57) 


Equations  (38)  and  (57)  can  be  substituted  into  equation  (54)  to 


obtain 


U. (1  -  %M  siny  ) 

sin8 - - - =— ± - 

c  uT(l  -  ^  siny^)  +  l^siny.^ 

Finally,  since  the  wavelets  and  packets  have  the  same  Initial 


direction  it  follows  that  y^  *  9c*  Thus  equation  (58)  can  be 
written  as 


M4sin46  -  7M3sin30  +  (18M2  -  8)sin20  -  20M  sin0  +8=0 

c  c  c  c 

(59) 

This  is  a  quartic  equation.  In  general,  the  expressions  for  the 
roots  of  such  an  equation  are  useless  for  computations  (Sokolnikoff 
and  Sokolnikoff,  1941).  One  thing  which  can  be  obtained  from -this 
equation  is  that  the  critical  angles  for  packets  are  a  function  of 
U^/G^  and  are  independent  of  period  and  shear  as  were  the  critical 
angles  for  the  wavelets. 

Often  equations  like  (54)  can  be  solved  for  specific  values 
by  iteration.  The  iteration  process  failed  in  this  case  due  to  non¬ 
convergence.  Finally  critical  angles  for  packets  were  determined  by 
trial  and  error  using  equations  (47)  and  (54) . 

Figures  2  and  3  show  graphically  the  critical  angles  obtained 
for  wavelets  and  packets.  From  these  figures  it  is  seen  that  packets 
reflect  at  smaller  angles  of  incidence  than  do  the  wavelets.  Thus 
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along  a  ray  the  packets  reflect  before  the  wavelets,  and  since  the 
rays  must  reflect  when  the  packets  reflect  (Section  II. D.)  it  follows 
that  the  ray  also  reflects  before  the  wavelets.  It  can  also  be  seen 
that  as  the  ratio  u^/G^  decreases  the  initial  angle  required  for 
a  ray  to  reflect  increases  and  eventually  the  critical  angles  for 
wavelets,  packets  and  rays  become  equal.  Further,  it  is  seen  that 
very  large  values  of  are  required  to  reflect  rays  of  small 

initial  angle.  For  wave  periods  consistent  with  surface  gravity 
waves  (7  to  20  seconds)  only  rays  with  initial  angles  greater  than 
44  degrees  will  be  reflected  by  open  ocean  currents.  To  reflect 
rays  with  initial  angles  less  than  44  degrees,  current  speeds  greater 
than  two  meters  per  second  must  be  encountered.  (The  maximum  speed 
of  the  Gulf  Stream  is  2  m/sec.) 

Figure  3  is  a  blow-up  of  a  section  of  Figure  2.  This  figure 
Indicates  the  range  of  critical  angles  most  likely  to  occur  for  sur¬ 
face  gravity  waves  in  the.  open  ocean. 

G.  Vorticity  equation 

Another  approach  to  determining  ray  paths  for  wave  groups  is 
to  employ  a  curvature  equation  similar  to  the  equation  discussed  by 
Kenyon  (1971)  and  used  by  Teague  (1974).  This  curvature  equation 
involves  the  current  vorticity  and  the  conventional  group  speed.  In 
this  section,  a  curvature  equation  for  rays  of  wave  packets  using 
geometric  group  speed  will  be  developed  and  compared  to  the  ray 
curvature  equation  mentioned  by  Kenyon  and  used  by  Teague.  These 
curvature  equations  will  be  compared  to  the  advection  equation  (30). 
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A  ray  curvature  equation  can  be  derived  from  equations  (6) 
and  (10)  assuming  that  u  <<  G'.  The  time  derivative  of  equation  (10) 
is 

P  «  G*(m/m)  +  u  (60) 


where  P  is  the  velocity  of  the  ray.  Multiplying  equation  (60)  by 
m  and  taking  the  derivative  of  the  result  with  respect  to  time  yields 


Equation  (.6)  can  be  expanded  and  rewritten  as 


dm/dt 


-  m  •  Vu  -  m  x  (V  x  u) 


just  as  was  done  with  equation  CIS)  for  the  wavelets. 

For  steady  state  conditions  3G'/9t  =  0  and  3u/3t  «=  0. 

•f 

Therefore,  along  a  ray  the  total  derivatives  for  G'  and  u  are 

dG'/dt  =*  P  •  VG*  C63) 


du/dt  =  P  •  Vu 


C64) 


When  (60)  is  substituted  into  (.64),  and  only  the  first  order  term  in 
u  is  retained,  it  is  found  that 

du/dt  *  G'(m/m)  •  Vu  (65) 


The  substitution  of  equations  (62),  (63),  and  (65)  into  equation  (61) 
yields 

d(Pm)/dt  -  -G’mVG*  -  G’m*  (V x  u)  +  m(?  •  VG')  +  u  dm/dt  (66) 
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Since  u  «  G'  it  can  be  assumed  that  P  -  G'  and  that  the  direction 
of  m  is  approximately  the  same  as  the  direction  of  the  ray.  Then 
P  can  be  substituted  for  G1  in  equation  (66)  and  a  unit  vector  in 
the  direction  of  the  ray,  defined  as  e^  •  P/P,  can  be  substituted  for 
the  direction  of  m.  Then 


d(Pm)/dt  =  S  d(Pm)/dt  +  Pm(der/dt)  (68) 

where  e^  and  de^/dt  are  perpendicular.  Equations  (67)  and  (68) 
can  be  combined  to  obtain 

d@^/dt  *  -VG'  -  e^x  (VxJi)  +  (u/Pm)  (dm/dt)  -  (er/Pm)d(Pm) /dt 

+  (er/P)(P  •  VG')  (69) 

The  third  term  on  the  right-hand  side  of  this  equation  can  be 
neglected  since  u  «  P  and  since  the  wave  properties  are  assumed 
to  change  little  over  distances  of  the  order  of  a  wavelength. 

The  direction  of  the  ray,  i.e.,  the  direction  of  e^  with 
respect  to  the  positive  x-axis,  is  denoted  by  p.  Then 

dp/dt  ■  nr  •  (de^/dt)  *  -n^  •  VG'  -  n^  •  £@rx(Vxu)J  (70) 

where  n^  is  a  unit  vector  perpendicular  to  e^. 

For  surface  gravity  waves,  only  the  vertical  component  C 
of  the  vorticity  affects  the  ray  trajectory.  The  z-component  of 
(7  x  u)  is  given  by 


vvs 


31 


C  *  3u  /3x  -  3u  /3y 

y  x 


Simplification  of  (70)  leads  to 


dp/dt  =  sinp  (3G'/3x)  -  cosp  (3G’/3y)  +  3u  /3x  -  3u  /3y 

y  x  (?: 

Further,  since  ds^  -  G’dt,  where  ds^  is  an  element  of  arc  length 
along  the  ray,  it  is  found  that 


l  I  ’ 

s  “  G* 
r  L 


1  r  <  „  3g'  „ 3o'  —i  3u* 

sinp  ^  -  cosp  +  3J*  - 


Equation  (73)  is  the  curvature  equation  for  computing  the  rays  of  wave 
packets  refracted  by  a  current.  This  equation  should  be  used  in  con¬ 
junction  with  the  wavelet  and  packet  curvature  equations. 

If  6  and  y  are  assumed  to  be  in  the  same  direction,  and 
if  the  group  speed  U*  is  assumed  to  be  constant  along  the  ray, 
equation  (73)  will  reduce  to  the  curvature  equation  discussed  by 
Kenyon  (1971)  and  used  by  Teague  (1974).  That  is 


,  f  3u  3u 
^  _ £  _ x 

U’  L3x  “  3y 


This  equation  will  determine  ray  paths,  but  only  under  very  strict 
limitations. 

The  group  speed  U’  is  not  constant  even  when  u  <<  U', 
since  a  current  causes  a  Doppler  shift  in  the  wave’s  frequency. 
However,  if  u/U’  is  very  small,  the  changes  in  U'  or  G'  along  a 
ray  will  be  quite  small.  Secondly,  0  and  y  are  not  generally  in 
the  same  direction  due  to  refraction,  but  their  difference  would  be 
small  when  u/U’  is  small. 
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A  Snell's  law  type  expression  can  be  determined  from  equation 
(73),  since  Snell's  law  is  the  integrated  form  of  the  ray  curvature 
equation.  To  integrate  equation  (73)  it  will  be  assumed  that  G'  is 
constant  along  a  ray,  that  the  current  has  parallel  speed  contours, 
and  that  a  Cartesian  coordinate  system  with  the  y-axis  parallel  to 
the  contours  is  used.  Thus  equation  (73)  can  be  expressed  as 

dp/dsr  ■  (1/G')  (du/dx)  (75) 

Separating  variables  and  substituting  cosp  for  dx/dsr  equation  (75) 
simplifies  to 

cosp  dp  =  du/G'  (76) 

Integrating  both  sides  of  this  equation  yields 

sinp  -  u/G'  **  constant  (77) 

This  result  Is  Snell's  law  for  rays,  and  It  can  be  expressed  as 

sinpi  -  u^/G'  =  sinpt  -  ut/G'  (78) 

If  u^  ■  0,  then 

sinpt  -  sinp^  +  u^G*  (79) 

When  turning  points  and  critical  angles  are  considered,  it 
is  found  that  the  critical  angle  equation  for  the  rays  has  the  form 

sinpc  -  1  -  M  (80) 

where  p  is  the  critical  angle  for  the  rays.  This  shows  that  the 


critical  angle  for  rays  is  a  function  of  u  /G  just  as  was  the 
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critical  angle  for  wavelets  and  packets. 

Since  critical  angles  for  rays  are  the  same  as  critical 
angles  for  packets,  it  should  be  possible  to  derive  equation  (.80) 
from  equation  (53)  when  u  <<  G'  and  when  G'  is  constant.  Thus, 
when  0  *  90  degrees,  G*  *  If'siny^  and  equation  (53)  can  be  rewritten 
as 

sin0  -  G' / (u  +  G')  C81) 

C  T 

and 

sin6  -  (1  +  u  /G')_1  (82) 

C  T 

Equation  (82)  can  be  expanded  to  obtain 

sin8  *  1  -  u  /G'  +  ,  ,  ,  higher  order  terms  (83) 

Thus,  neglecting  the  higher  order  terms,  equation  (53)  reduces  to 
equation  (80). 

Physical  properties  in  addition  to  those  discussed  in 
Section  III.D.  become  evident  when  equation  (73)  is  viewed  in  the 
form 

dp/dsf  -  5/G’  (84) 

From  this  equation  it  can  be  seen  that  the  curvature  of  a  ray  has  the 
same  sign  as  the  vorticity  of  the  current.  Therefore,  for  a  variable 
following  (opposing)  current,  rays  bend  toward  decreasing  (increasing) 
current  speed.  Further,  under  the  right  conditions  rays  in  a  follow¬ 
ing  (opposing)  current  may  be  trapped  about  a  local  minimum  (maximum) 


In  current  speed.  It  can  also  be  seen  that  the  curvature  of  a  ray 
decreases  with  increasing  geometric  group  speed  (i.e.,  Increasing 
period).  So  small  period  waves  are  affected  to  a  greater  extent  by 
a  current  than  are  large  period  waves.  These  properties  appear  to 
occur  (.Section  IV. C.)  even  when  the  assumptions  made  in  obtaining 
equation  (84)  do  not  hold. 
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IV.  CURRENT  MODELS  AND  REFRACTION  RESULTS 


A.  Description  of  current  models 

The  refraction  of  wave  packets  by  currents  will  be  examined 
for  four  current  models.  Two  of  the  current  models  will  be  patterned 
after  the  Gulf  Stream.  The  other  two  models  will  be  patterned  after 
the  Circumpolar  Current.  The  Gulf  Stream  was  chosen  since  it  is  the 
swiftest  deep  ocean  current  and  since  its  effect  on  waves  could 
greatly  influence  wave  patterns  off  the  east  coast  of  the  United 
States.  The  Circumpolar  Current  was  chosen,  since  for  the  waves 
being  considered  its  current  speed  to  geometric  group  speed  ratio  is 
relatively  small  (<  0.1),  and  since  it  affects  a  large  portion  of  the 
Atlantic  and  Pacific  Oceans. 

Two  parallel  current  models  are  considered.  They  are  assumed 
to  be  unidirectional  shear  currents,  uniform  in  the  direction  of  flow, 
and  Independent  of  depth.  One  of  the  parallel  current  models  has  a 
maximum  current  speed  of  50  cm/sec  and  a  width  of  1000  km.  These 
values  are  representative  of  the  Circumpolar  Current.  The  current 
speed  is  assumed  to  have  a  parabolic  profile  such  that  its  speed 
Increases  from  zero  at  the  edge  of  the  current  to  its  maximum  value 
at  the  center  of  the  current.  The  speed  variation  through  the 
current  is  expressed  by 

u  -  -  2  x  10_4(W2  -  1000  W)  (85) 


ilfli'A'. 


k’  !%m'.  V.  /U'i  ,~W«.  . 
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where  u  is  the  speed  of  the  current  in  cm/sec,  and  W  is  the 
cp 

distance  from  the  edge  of  the  current  in  km.  The  results  obtained 
from  this  equation  are  plotted  in  Figure  4. 

The  parallel  Gulf  Stream  model  has  a  maximum  current  speed 
of  2  m/sec  and  a  width  of  100  km.  This  current  model  also  has  a 
parabolic  profile  such  that  the  current  speed  increases  from  zero 
at  the  edge  of  the  current  to  its  maximum  value  at  the  center  of  the 
current.  The  speed  across  this  current  is  expressed  by 

u  -  -  8  x  10_2(W2  -  100  W)  (.86) 

gs 

where  u  is  the  speed  of  the  current  in  cm/sec.  The  results  obtained 
gs 

from  this  equation  are  shown  in  Figure  5. 

The  curved  current  models  have  the  same  current  speed  profiles 
and  dimensions  as  the  parallel  currents,  and  are  assumed  to  be  inde¬ 
pendent  of  depth.  However,  their  contours  are  assumed  to  be  sections 
of  horizontal,  non-divergent  circular  rings.  The  Circumpolar  Current 
is  constructed  so  that  its  inner  radius  of  curvature  is  4000  km  and 
its  outer  radius  of  curvature  is  5000  km,  similar  to  Kenyon's  model. 

The  curved  Gulf  Stream  model  has  an  inner  radius  of  curvature  of 
500  km  and  an  outer  radius  of  curvature  of  600  km,  which  roughly 
corresponds  to  the  Gulf  Stream  off  the  coast  of  Georgia.  To  examine 
the  refraction  effect  of  non-parallel  current  contours  the  wave 
trajectories  obtained  for  the  curved  current  models  can  be  compared 
with  the  rays  determined  for  the  corresponding  parallel  current 
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B.  Wave  period 

In  this  study,  waves  of  periods  from  seven  to  twenty  seconds 
are  considered.  This  is  the  significant  range  of  periods  found 
coastward  of  the  Gulf  Stream.  Figures  6,  7,  and  8  show  the  annual 
significant  wave  period  distribution  measured  at  sites  off  the  east 
coast  of  the  United  States  (Thompson,  1977).  It  is  seen  from  these 
figures  that  swell  wave  periods  peak  near  10  seconds.  In  the  Southern 
Ocean  there  are  frequently  large  storms  and  a  large  fetch.  As  a 
result,  longer  period  swell  waves  (.15  to  20  seconds)  are  expected  to 
be  more  prevalent  in  the  Circumpolar  Current  than  in  the  Gulf  Stream. 

C.  Properties  of  ray  trajectories: 
parallel  and  curved  contours 

In  this  section  ray  trajectories  computed  from  equations  (.26) , 
(28),  and  (30),  are  examined  for  waves  in  the  current  models  previously 
described.  For  each  model  and  wave  period  five  Initial  angles  are 
considered.  The  initial  angles  are  15,  30,  45,  60,  and  75  degrees 
measured  relative  to  the  normal  to  the  current  contours. 

The  ray  properties  deduced  from  the  theory  in  Section  III.D. 
can  be  demonstrated  through  these  trajectories.  It  was  determined 
that  p  >  8  for  a  following  current,  that  p  <  9  for  an  opposing 
current,  and  that  p  ■  0  ■  90  degrees  at  a  turning  point.  These 
properties  are  presented  graphically  in  Figures  9  thru  13  which  show 
the  variation  in  p,  8,  and  y  along  some  representative  rays. 

Figure  9  depicts  the  values  of  the  angles  along  a  ray  in  a  following 
current,  and  Figure  10  is  for  a  ray  in  an  opposing  current.  It  is 


Figure  7,  Annual  significant  period  distribution  from  Holden  Beach, 
North  Carolina;  computed  from  1,024-second  digital  wave 
records  taken  four  times  daily.  (Thompson,  1977) 


Figure  8.  Annual  significant  period  distributions  from  Wrightsville 
Beach,  North  Carolina;  computed  from  1,024-second  digital 
wave  records  taken  four  times  daily.  (Thompson,  1977) 
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clearly  seen  that  for  the  following  current  p  >  9,  and  that  p  <  0 
for  the  opposing  current.  It  is  also  seen  that  the  angle  y  changes 
the  least  of  the  three  angles.  This  is  because  the  speed  of  the  wave¬ 
lets  is  greater  than  that  of  the  packets  and  rays.  Thus,  the  wavelets 
are  affected  the  least  hy  the  current.  Figure  11  is  an  example  of 
the  variation  in  the  angles  along  a  ray  which  goes  through  a  turning 
point.  It  is  seen  that  when  a  ray  totally  reflects  the  values  of  p 
and  9  steadily  increase,  and  reach  90  degrees  at  the  turning  point. 
The  angle  y,  however,  is  seen  to  increase  to  its  maximum  value  at 
the  turning  point  and  then  to  decrease  to  its  initial  value.  Due  to 
the  symmetry  of  the  current  models  used  p,  0,  and  y  have  the  same 
value  when  the  ray  exits  the  current  as  when  it  entered  the  current, 
except  for  p  and  8  when  the  ray  totally  reflects. 

It  was  also  determined  that  when  u  «  G*  the  values  of  p 

and  0  are  nearly  equal.  In  fact,  it  was  found  that  in  such  cases 

y  is  also  nearly  equal  to  p  and  9.  An  example  of  this  is  seen 

in  Figure  12  where  N  -  u  /G.  ■  0.046.  For  this  ray  p  and  0 

max  i 

differ  by  no  more  than  1.5  degrees  and  0  and  Y  differ  by  no  more 
than  1  degree.  However,  in  Figure  9,  where  N  *  0.366,  p  and  0 
differ  by  about  12  degrees,  and  0  and  y  differ  by  about  4  degrees. 
This  changes  somewhat  if  the  ray  totally  reflects  (Figure  13) .  When 
there  is  a  reflection  p  and  0  remain  nearly  equal,  but  y  deviates 
considerably  from  the  other  angles. 

Further,  at  a  turning  point,  it  was  determined  that  either 
u  ■  G*  ■  0  or  u  ■  -G’ .  The  first  case  is  seen  to  occur  for  rays 

X  X  X  X 

which  totally  reflect  in  a  following,  parallel  current.  A  good 


example  of  this  is  presented  in  ray  5,  Figure  14.  At  the  turning 
point  the  ray,  packet,  and  current  directions  become  parallel.  The 
second  case  occurs  for  curved  currents.  It  was  found  that  the 
x-component  of  current  and  wave  packet  velocities  balanced  each 
other  in  rays  4  and  5,  Figures  26  and  30.  Also,  the  current  and 
ray  directions  were  found  to  directly  oppose  each  other.  An 
interesting  observation  is  that  these  rays  reflect  in  an  opposing 
current  rather  than  in  a  following  current. 

Lastly,  it  was  deduced  in  Section  III.D.  that  in  an  opposing 
current  rays  can  be  turned  so  they  would  eventually  become  perpen¬ 
dicular  to  the  current  contours.  This  is  seen  in  Figures  19  and  21, 
which  show  rays  in  the  Gulf  Stream.  However,  for  the  current  model 
with  the  slower  current  speeds  C 50  cm/sec  maximum)  there  is  not 
sufficient  refraction  to  cause  the  rays  to  become  perpendicular  to 
the  contours  (Figures  15  and  17). 

Some  ray  properties  were  deduced  in  Section  III.G.  from  the 
vortlclty  equation  which  is  an  approximate  equation  useful  only 
when  u  «  G'.  Even  so,  it  is  apparent  from  the  trajectories  that 
these  properties  are  general  and  occur  even  when  u/G'  is  as  large 
as  0.366.  In  accordance  with  equation  (84)  rays  are  bent  away  from 
the  current  normal  in  a  following  current,  and  toward  the  current 
normal  in  an  opposing  current.  The  rays  in  Figures  18  and  19  clearly 
demonstrate  these  results.  Also,  it  was  determined  that  the  curva¬ 
ture  of  a  ray  increases  as  the  vortlclty  increases  and  as  the  wave 
period  decreases.  The  average  vortlclty  of  the  current  shown  in 
Figure  22  is  1  x  10  ®  sec”*.  The  average  vorticity  of  the  current 
shown  in  Figure  18  is  4  x  10”^  sec  *. 


The  wave  period  is  seven 


seconds  for  the  rays  in  both  figures.  It  is  clear  that  the  curva¬ 
ture  of  the  rays  shown  in  Figure  18  is  greater  than  the  curvature 
of  the  rays  shown  in  Figure  22.  In  Figures  18  and  20  the  average 
vorticity  of  the  current  is  4  *  10  ^  sec  The  wave  periods, 
however,  differ  in  these  figures.  Seven  second  period  wave  rays 
are  depicted  in  Figure  18  and  10  second  period  wave  rays  are  shown 
in  Figure  20.  The  rays  with  the  smaller  period  are  seen  to  have  the 
greater  curvature. 

Ray  trajectories  computed  for  the  parallel  current  models 
are  shown  in  Figures  14  thru  21.  For  the  parallel  current  patterned 
after  the  Circumpolar  Current,  wave  periods  of  14  and  17  seconds  were 
considered.  These  periods  correspond  to  N  values  of  0.046  and 
0.038,  respectively.  In  Figures  14  thru  17  it  is  seen  that  little 
refraction  occurs  for  these  values  of  N.  It  is  also  seen  that  only 
rays  with  large  Initial  angles  are  totally  reflected,  as  was  indi¬ 
cated  in  Figure  2.  The  critical  angle  computed  when  N  is  0.046  is 
71  degrees.  When  N  is  0.038  the  critical  angle  is  74  degrees.  The 
rays  shown  in  Figures  14  and  16  are  consistent  with  these  results. 

Ray  trajectories  for  the  current  model  patterned  after  the 
Gulf  Stream  are  shown  in  Figures  18  thru  21.  Wave  periods  of  7  and 
10  seconds  are  considered.  These  periods  correspond  to  N  values  of 
0.366  and  0.256,  respectively.  The  rays  in  this  current  are  seen  to 
undergo  considerable  refraction.  The  critical  angles  computed  for 
these  rays  are  44  degrees  when  N  is  0.366  and  50  degrees  when  N 
is  0.256. 

Ray  trajectories  computed  using  the  curved  current  models  are 
shown  in  Figures  23  thru  38.  It  is  evident  from  these  figures  that 


the  trajectories  for  waves  in  curved  currents  can  be  considerably 
different  from  the  trajectories  for  waves  in  parallel  currents  even 
though  the  current  speed  profile  and  current  dimensions  are  the 
same.  This  can  be  explained  by  the  greater  variation  in  current 
vorticity  which  occurs  for  curved  currents  compared  to  parallel 
currents. 

With  curved  currents  the  side  from  which  the  ray  enters  the 
current  becomes  important.  With  parallel  currents  it  is  only  necessary 
to  distinguish  between  following  and  opposing  currents,  since  the 
symmetry  of  the  currents  results  in  the  same  trajectories  no  matter 
from  which  side  of  the  current  the  waves  enter.  With  the  curved 
currents  this  is  not  so.  The  trajectoreis  of  waves  which  enter  from 
the  inner  radius  side  of  the  current  are  different  from  the  trajectories 
which  enter  from  the  outer  radius  side.  This  can  he  seen  in  Figures 
23  and  24  and  in  Figures  31  and  32.  It  is  also  seen  that  the  critical 
angles  are  different  for  the  inner  radius  side  than  for  the  outer 
radius  side.  For  example,  in  Figure  23  the  critical  angles  for  these 
rays  must  be  greater  than  75  degrees.  Zn  Figure  24,  the  critical 
angle  is  found  to  be  less  than  60  degrees.  A  similar  result  is  seen 
in  the  other  figures.  Thus,  the  critical  angle  for  waves  which  enter 
a  curved  current  from  the  inner  radius  side  is  larger  than  the  critical 
angle  for  waves  which  enter  from  the  outer  radius  side  given  the  same 
wave  period  and  current  speed  profile.  These  critical  angles  differ 
from  those  computed  for  parallel  currents  (Figure  2) .  The  critical 
angles  for  rays  which  enter  from  the  inner  radius  side  are  greater 
than  the  critical  angles  determined  for  the  same  waves  in  a  parallel 
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current.  The  critical  angles  for  rays  which  enter  from  the  outer 
radius  side  are  less  than  those  determined  for  parallel  currents. 

The  refraction  of  wave  packets  by  curved  currents  can 
result  in  some  very  interesting  rays.  For  example,  rays  may,  due  to 
the  configuration  of  the  current,  pass  through  the  current  and 
then  re-enter  the  current  further  downstream  or  upstream.  This  is 
seen  to  occur  for  ray  3  in  Figures  26  and  30.  This  can  occur  for 
either  a  following  or  an  opposing  current,  but  only  if  the  ray  exits 
the  current  on  the  inner  radius  side.  It  can  also  be  observed  from 
these  figures  that  rays  may  go  through  several  inflection  points 
Cray  4) ,  and  that  rays  may  be  focused  into  a  particular  area  (rays  4 
and  5) .  Comparison  of  Figures  26  and  30  indicates  that  the  area  into 
which  the  rays  are  focused  moves  further  away  from  the  entry  point 
of  the  waves  into  the  current  as  the  wave  period  of  the  rays 
Increases.  One  might  say  that  the  focal  length  increases  as  the 
wave  period  Increases.  .  Very  large  period  waves  (.-  20  seconds)  would 
be  focused  outside  of  the  current. 

Ray  3,  Figure  32;  rays  4  and  5,  Figures  34  and  38;  and  ray  4, 
Figure  25  are  seen  to  abruptly  stop  within  the  current.  This  occurs 
due  to  numerical  problems  which  come  about  when  G'  goes  to  zero. 
After  a  total  reflection,  it  has  been  shown  that  3  and  y  differ 
considerably.  In  these  rays,  this  difference  results  in  G '  becoming 
zero.  Further  study  will  he  needed  to  determine  what  happens  to  the 
ray  when  this  occurs. 
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Figure  9 


The  variation  in  the  values  of  the  angles  p,  9,  and  Y 
along  ray  2  in  Figure  18.  This  ray  has  an  initial  angle 
of  30  degrees,  and  passes  through  a  following  current 
without  reflection.  (N  -  0.366) 
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Figure  11.  The  variation  in  the  values  of  the  angles  p,  8,  and  y 
along  ray  4  in  Figure  18.  This  ray  has  an  initial  angle 
of  60  degrees.  This  ray  is  reflected.  (N  -  0.366) 
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Figure  13.  The  variation  in  the  values  of  the  angles  p,  6,  and  y 
along  ray  5  in  Figure  14.  This  ray  has  an  initial  angle 
of  75  degrees  and  is  totally  reflected.  (N  -  0.046) 
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Figure  14.  Ray  trajectories  for  14  second  period  waves  in  a 
following  current.  The  current  model  used  was 
the  parallel  model  patterned  after  the  Circum¬ 
polar  Current.  Current  speeds  are  shown  in  cm/sec 
N  ■  0.046,  and  the  positive  y-axis  is  in  the 
direction  of  the  current.  Grid  spacing  is  20  km. 
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Figure  15.  Ray  trajectories  for  14  second  period  waves  In  an 
opposing  current.  The  current  model  used  was  the 
parallel  model  patterned  after  the  Circumpolar 
Current.  Current  speeds  are  shown  In  cm/sec, 

N  «  0.046,  and  the  positive  y-axis  is  in  the  direc 
tion  of  the  current.  Grid  spacing  is  20  km. 
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Figure  16.  Ray  trajectories  for  17  second  period  waves  in  a 

following  current.  The  current  model  used  was  the 
parallel  model  patterned  after  the  Circumpolar 
Current.  Current  speeds  are  shown  in  cm/sec, 

N  ■  0.038,  and  the  positive  y-axis  is  in  the  direc¬ 
tion  of  the  current.  Grid  spacing  is  20  km. 


Figure  17,  Ray  trajectories  for  17  second  period  waves  in  an 
opposing  current.  The  current  model  used  was  the 
parallel  model  patterned  after  the  Circumpolar 
Current.  Current  speeds  are  shown  in  cm/sec, 

N  ■  0.038,  and  the  positive  y-axis  is  in  the  direc¬ 
tion  of  the  current.  Grid  spacing  is  20  km. 


Figure  19.  Ray  trajectories  for  7  second  period  waves  in  an 
opposing  current.  The  current  model  used  was  the 
parallel  model  patterned  after  the  Gulf  Stream. 
Current  speeds  are  shown  in  cm/sec,  N  -  0.366, 
and  the  positive  y-axis  is  in  the  direction  of  the 
current.  Grid  spacing  is  2.5  km. 


50 


200 


50 


Figure  21.  Ray  trajectories  for  10  second  period  waves  in  an 
opposing  current.  The  current  model  used  was  the 
parallel  model  patterned  after  the  Gulf  Stream. 
Current  speeds  are  shown  in  cm/sec,  N  -  0.256, 
and  the  positive  y^axis  is  in  the  direction  of 
the  current.  Grid  spacing  is  2.5  km. 


Figure  22.  Ray  trajectories  for  7  second  period  waves  in  a 

following  current.  The  current  model  used  was  the 
parallel  model  patterned  after  the  Circumpolar 
Current.  Current  speeds  are  shown  in  cm/sec, 

N  ■  0.091,  and  the  positive  y-axis  is  in  the  direc 
tion  of  the  current.  Grid  spacing  is  20  km. 


Figure  23.  Ray  trajectories  for  14  second  period  waves  in  the  Cir¬ 
cumpolar  Current  model.  The  waves  enter  a  following 
current  from  the  inner  radius  side.  Current  speeds  are 
given  in  cm/ sec  and  N  *  0.046.  Grid  spacing  is  20  km. 


Figure  24.  Ray  trajectories  for  14  second  period  waves  in  the  Circum¬ 
polar  Current  model.  The  waves  enter  a  following  current 
from  the  outer  radius  side.  Current  speeds  are  given  in 
cm/sec  and  N  -  0.046.  Grid  spacing  is  20  km. 


Figure  25.  Ray  trajectories  for  14  second  period  waves  in  the  Cir¬ 
cumpolar  Current  model.  The  waves  enter  an  opposing 
current  from  the  inner  radius  side.  Current  speeds  are 
given  in  cm/sec  and  N  «  0.046.  Grid  spacing  is  20  km. 


Figure  26.  Ray  trajectoreis  for  14  second  period  waves  in  the  Cir¬ 
cumpolar  Current  model.  The  waves  enter  an  opposing 
current  from  the  outer  radius  side.  Current  speeds  are 
given  in  cm/ sec  and  N  *  0.046.  Grid  spacing  is  45  km. 


Figure  28.  Ray  trajectories  for  17  second  period  waves  in  the  Cir¬ 
cumpolar  Current  model.  The  waves  enter  a  following 
current  from  the  outer  radius  side.  Current  speeds  are 
given  in  cm/sec  and  N  -  0.038.  Grid  spacing  is  20  km. 


Figure  29.  Ray  trajectories  for  17  second  period  waves  in  the  Cir¬ 
cumpolar  Current  model.  The  waves  enter  an  opposing 
current  from  the  inner  radius  side.  Current  speeds  are 
given  in  cm/sec  and  N  •  0.038.  Grid  spacing  is  20  km. 
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Figure  30.  Ray  trajectories  for  17  second  period  waves  in  the  Cir¬ 
cumpolar  model.  The  waves  enter  an  opposing  current 
from  the  outer  radius  side.  Current  speeds  are  given 
in  cm/sec  and  N  ■  0.038.  Grid  spacing  is  45  km. 


Figure  31.  Ray  trajectories  for  7  second  period  waves  in  the  curved 
Gulf  Stream  model.  The  waves  enter  a  following  current 
from  the  inner  radius  side.  Current  speeds  are  given  in 
cm/sec  and  N  ■  0.366.  Grid  spacing  is  2.5  km. 


Figure  32.  Ray  trajectories  for  7  second  period  waves  in  the  curved 
Gulf  Stream  model.  The  waves  enter  a  following  current 
from  the  outer  radius  side.  Current  speeds  are  given  in 
cm/sec  and  N  -  0.366.  Grid  spacing  is  2.5  km. 


Figure  33.  Ray  trajectories  for  7  second  period  waves  in  the 

curved  Gulf  Stream  model.  The  waves  enter  an  opposing 
current  from  the  inner  radius  side.  Current  speeds 
are  in  cm/sec  and  N  ■  0.366.  Grid  spacing  is  2.5  km. 


Figure  35.  Ray  trajectories  for  10  second  period  waves  in  the  curved 
Gulf  Stream  model.  The  waves  enter  a  following  current 
from  the  inner  radius  side.  Current  speeds  are  in  cm/sec 
and  N  -  0.256.  Grid  spacing  is  2.5  km. 
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Figure  37.  Ray  trajectories  for  10  second  period  waves  in  the  curved 
Gulf  Stream  model.  The  waves  enter  an  opposing  current 
from  the  inner  radius  side.  Current  speeds  are  given  in 
cm/sec  and  N  -  0.256.  Grid  spacing  is  2.5  km. 
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Figure  38.  Ray  trajectories  for  10  second  period  waves  in  the 

curved  Gulf  Stream  model.  The  waves  enter  an  opposing 
current  from  the  outer  radius  side.  Current  speeds  are 
given  in  cm/sec  and  N  *  0.256.  Grid  spacing  is  2.5  km 


V.  COMPARISON  OF  REFRACTION  METHODS 

A.  Monochromatic  refraction 

Arthur  (.1950)  considered  the  refraction  of  monochromatic 
waves  by  currents  in  shallow  water.  The  equations  he  used  and  the 
results  he  obtained  were  discussed  in  Chapter  II.  In  this  section 
Arthur's  approach  is  applied  to  deep  water  waves. 

Monochromatic  trajectories  were  computed  using  equation  (26) 
which  determines  the  direction  of  the  waves  along  a  ray,  and 

tanp  ’  V'cosY  +  tanY  (87) 

which  determines  the  direction  of  the  ray.  These  equations  are  the 
same  as  equations  (.2)  and  (3)  adjusted  for  deep  water.  Comparison 
of  equations  (33),  for  wave  packets,  and  (87),  for  monochromatic 
waves,  indicates  that  the  ray  directions  for  monochromatic  waves 
changes  less  than  the  ray  direction  for  wave  packets,  since  G*  <  V' 
Also,  from  these  equations  it  is  found  that  the  monochromatic  rays 
and  monochromatic  waves  totally  reflect  at  the  same  point,  and  that 
their  critical  angles  are  larger  than  the  critical  angles  for  packet 
rays  (same  as  in  Figure  2).  These  characteristics  are  clearly  seen 
in  Figures  39  thru  41,  where  the  subscript  m  denotes  monochromatic 
rays  and  the  subscript  p  denotes  packet  rays.  Since  the  ray 
equation  for  monochromatic  rays  has  the  same  form  as  does  the  ray 
equation  for  wave  packets  the  same  general  physical  properties  apply 
to  both. 
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In  Sections  III.D.  and  IV, C.,  it  was  shown  that  p,  9,  and 
Y  are  nearly  equal  when  u  <<  G',  since  then  little  refraction 
occurs.  This  means  that  as  u/G'  approaches  zero  wave  packet  rays 
and  monochromatic  rays  become  nearly  identical.  However,  they  do 
not  become  identical,  since  G*  is  not  equal  to  V*.  The  best  agree¬ 
ment  between  monochromatic  rays  and  wave  packet  rays  is  seen  to 
occur  for  the  rays  with  small  initial  angles,  with  differences 
increasing  as  the  initial  angle  Increases.  The  worst  agreement  is 
for  rays  which  undergo  total  reflection.  This  occurs  since  the 
difference  in  the  angles  6  and  y  increases  as  the  initial  angle 
increases. 

B.  Vorticity  approach 

The  vorticity  equation  C74)  developed  in  Section  III.G. 
is  the  same  equation  used  By  Teague  (1974)  to  compute  ray  paths  through 
eddy  type  currents.  This  equation  was  derived  by  asstiming  that 
u  «  G' ,  that  U*  is  constant  along  a  ray,  and  that  the  packet  and 
wavelet  directions  are  identical.  Although  U'  is  not  actually  con¬ 
stant  along  a  ray,  it  changes  so  slightly  when  u  «  G'  that  there 
is  no  deceroable  difference  between  rays  computed  with  U'  constant 
and  rays  with  U*  allowed  to  vary.  It  was  found  that  when  N  was 
0.091,  the  conventional  group  speed  changed  by  less  than  4  percent 
along  a  ray  and  by  less  than  2  percent  when  N  was  0.046. 

Rays  computed  using  the  vorticity  equation  (74)  are  denoted 
in  the  figures  by  the  subscript  v.  It  is  clearly  seen  from  these 
figures  that  the  best  agreement  between  the  vorticity  rays  and  the 


wave  packet  rays  occurs  when  N  is  less  than  0.05  (Figure  41), 
and  when  the  initial  angles  are  small.  As  N  increases  the  vor- 
ticity  rays  become  less  reliable,  with  the  worst  agreement  in 
Figure  39  occurring  when  N  is  0.366.  Vorticity  rays  and  packet 
rays  are  seen  to  be  identical  or  nearly  so  when  the  initial  angles 
are  15  degrees  or  less,  and  to  be  significantly  different  when  the 
initial  angles  equal  or  exceed  the  critical  angles. 

The  vorticity  rays  shown  in  Figure  39  deviate  appreciably 
from  packet  rays,  since  the  assumptions  made  in  deriving  the  vorticity 
equation  (74)  break  down  for  such  high  values  of  N.  In  this  case 
the  speeds  U'  and  G*  ■  U'coaXO-y)  change  significantly  along  a  ray. 
It  was  found  that  when  N  «*  0.366,  U'  changes  by  about  5  percent  for 
a  ray  with  an  initial  angle  of  15  degrees  and  by  as  much  as  13  percent 
for  a  ray  with  an  initial  angle  of  45  degrees.  The  speed  G'  changes 
by  5  percent  along  a  ray  with,  an  initial  angle  of  15  degrees  to  as 
much  as  56  percent  along  a  ray  with  an  initial  angle  of  60  degrees. 

Teague  (1974)  considered  the  refraction  of  wave  packets  by 
currents  using  the  vorticity  method.  He  investigated  cyclonic 
eddies  and  swell  waves  with  N  values  of  0.17,  0.08,  and  0.02. 

From  what  has  been  shown  here,  it  is  evident  that  the  vorticity 
method  would  result  in  reasonable  ray  trajectories  for  the  second 
and  third  cases,  but  questionable  trajectories  for  the  first  case. 

C.  Kenyon's  approach 

Kenyon  (1971)  considered  the  refraction  of  wave  packets  by 
assuming  that  the  packets  travel  in  the  same  direction  as  the  wavelets 


within  the  packet.  Thus,  he  used  the  conventional  group  speed  instead 
of  the  geometric  group  speed.  Kenyon  determined  ray  trajectories  by 
integrating  the  ray  equations  (.4)  for  specific  currents.  He  stated 
that  his  approach  would  be  similar  to  the  vorticity  approach  when 
u  «  G'. 

A  method  of  determining  ray  trajectories,  which  yields  the 
same  results  as  was  obtained  by  Kenyon  (1971) ,  is  to  use  equation  (26) 
to  determine  the  wavelet  direction  along  a  ray  and 

tanp  “  u^i7  +  tanY  (88) 

to  compute  the  ray  direction.  Trajectories  computed  using  these 
equations  are  shown  in  Figures  39  thru  41.  They  are  denoted  by  the 
subscript  k.  Comparison  of  equation  (33)  and  equation  (88)  indicates 
that  rays  computed  using  these  equations  will  differ  as  0  and  y 
differ.  It  is  also  evident  from  these  equations  that  the  general 
physical  properties  discussed  in  Section  III.D.  and  III.G.  will  also 
apply  to  rays  computed  using  Kenyon's  method.  These  features  can  be 
seen  in  the  rays  shown  in  the  figures. 

It  is  apparent  from  the  figures  that  when  N  K  0.05,  Kenyon's 
method  and  the  one  derived  here  give  very  similar  rays.  But,  as  N 
Increases  the  difference  in  the  ray  trajectories  increases,  since 
then  the  difference  in  the  directions  of  the  wavelets  and  packets 
hecomes  significant.  A  similar  relationship  is  found  for  initial 
angles.  With  small  initial  angles  there  is  better  agreement  than 
for  larger  initial  angles. 

Lastly,  it  can  be  seen  that  when  u  «  G'  the  vorticity  method 
gives  rays  very  similar  to  those  determined  by  Kenyon's  method. 


D.  Summary  of  refraction  methods 

In  this  chapter,  it  has  been  shown  that  the  vorticity  method 
is  valid  only  under  very  limited  conditions,  with  the  best  results 
occurring  when  N  £  0.05.  Reasonable  trajectories,  though,  can  be 
obtained  even  if  N  is  as  large  as  0.10.  Thus,  it  can  be  concluded 
that  for  swell  waves  (7  to  20  seconds)  the  vorticity  method  of  wave 
refraction  would  not  be  applicable  to  currents  like  the  Gulf  Stream 
(N  ranges  from  0.366  to  0,13),  but  could  be  applied  to  currents  like 
the  Circumpolar  Current  (N  ranges  from  0.09  to  0.03).  The  vorticity 
method  should,  therefore,  be  confined  to  currents  with  speeds  less 
than  about  50  cm/sec  for  swell  waves. 

The  results  presented  in  this  chapter  indicate  that  the 
refraction  methods  discussed  yield  similar  ray  trajectories  when 
u  «  G’  and  when  the  initial  angles  are  small  (<  30  degrees).  It 
has  also  been  shown  that  as  N  increases  and  as  the  initial  angle 
Increases  the  ray  trajectories  become  quite  different  using  the 
various  refraction  methods.  It  can  be  seen  from  Figures  39  thru  41 
that  rays  computed  using  the  monochromatic  method,  vorticity  method, 
or  Kenyon's  method  agree  more  consistently  with  each  other  than  any 
of  them  do  with  rays  computed  using  the  wave  packet  method  presented 
in  this  study.  In  only  the  latter  method  is  the  difference  in  the 
wavelet  and  packet  directions  taken  into  account. 
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Figure  39.  The  rays  in  this  figure  were  computed  using  the 
various  refraction  methods.  The  subscript  p 
denotes  the  wave  packet  refraction  method  developed 
in  this  paper.  The  subscript  k  denotes  Kenyon's 
method  of  wave  packet  refraction  while  the  subscripts 
v  and  m  denote  the  vorticity  and  monochromatic  methods, 
respectively.  Rays  with  initial  angels  of  15,  30, 

45,  and  60  degrees  are  shown.  Current  speeds  are 
given  in  cm/sec  and  N  ■  0.366.  Grid  spacing  is  2.5  km. 
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Figure  40.  The  rays  in  this  figure  were  computed  using  the  various 
refraction  methods.  The  subscript  p  denotes  the  wave 
packet  refraction  method  developed  in  this  paper.  The 
subscript  k  denotes  Kenyon's  method  of  wave  packet 
refraction  while  the  subscripts  v  and  m  denote  the  vor- 
ticity  and  monochromatic  methods,  respectively.  Rays 
with  Initial  angles  of  15,  30,  45,  60,  and  75  degrees 
are  shown.  Rays  lv,  lk>  and  lp  are  seen  to  be  identical 
as  are  rays  5m  and  5^.  Current  speeds  are  given  in  cm/sec 
and  N  ■  0.091.  Grid  spacing  is  20  km. 


The  rays  in  this  figure  were  computed  using  the  various 
refraction  methods.  The  subscript  p  denotes  the  wave 
packet  refraction  method  developed  in  this  paper.  The 
subscript  k  denotes  Kenyon’s  method  of  wave  packet 
refraction  while  the  subscripts  v  and  m  denote  the  vor- 
ticity  and  monochromatic  methods,  respectively.  Rays 
with  initial  angles  of  15,  30,  45,  60,  and  75  degrees 
are  shown.  Rays  lv,  1^,  and  lp  are  identical  as  are 
rays  2V  and  2p.  Current  speeds  are  given  in  cm/sec 
and  N  ■  0.046.  Grid  spacing  is  20  km. 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  Conclusions 

Breeding's  method  (.1972,  1978,  and  1981)  of  wave  packet 
refraction  has  been  applied  in  this  study  to  the  refraction  of  wave 
packets  by  currents.  In  this  method,  the  direction  of  the  wavelets 
and  of  the  packets  vary  after  refraction.  Thus,  curvature  expressions 
for  both  wavelets  and  packets  are  needed  for  ray  trajectory  computa¬ 
tions.  It  was  found  from  the  expression  for  the  wavelet  direction 
along  a  ray  (equations  21  and  26)  that  the  change  in  direction  of 
the  wavelets  is  dependent  primarily  upon  the  shear  of  the  current  and 
the  phase  speed  relative  to  a  fixed  point.  The  change  in  direction 
of  the  packets  (equations  27  and  28)  was  found  to  be  dependent  upon 
the  current  shear,  the  geometric  group  speed  relative  to  a  fixed 
point,  the  direction  of  the  wavelets  y,  and  the  difference  in  the 
direction  between  the  wavelets  and  packets  $. 

Arthur  (1950)  and  Kenyon  (1971)  pointed  out  that  in  a  current 
rays  are  not  generally  orthogonal  to  the  wave  fronts.  Thus,  the 
direction  of  the  ray  is  different  from  the  direction  of  either  the 
packets  or  the  wavelets  along  the  ray.  It  was  shown  that  an  advectlon 
equation  (equation  30)  patterned  after  Arthur's  approach,  or  a  vorti- 
clty  equation  (equation  73  or  74)  could  be  derived  to  determine  ray 
paths.  The  advectlon  equation  was  found  to  be  more  general  and  more 
accurate  than  the  vortlcity  equation.  The  best  results  from  the 


vorticity  equation  can  be  expected  when  N  £  0.0S.  Values  of  N  as 
large  as  0.10  would  give  reasonable  results,  but  N  values  greater 
than  this  would  be  questionable.  Thus,  it  is  concluded  that  a  set 
of  three  equations  is  required  to  compute  ray  trajectories  for  wave 
packets  refracted  by  a  current.  The  ray  is  determined  by  simultaneously 
solving  the  curvature  expression  for  the  wavelets  (equation  26) .  the 
curvature  expression  for  the  packets  (equation  28) .  and  the  advection 
equation  for  the  rays  (equation  30) . 

The  curvature  equations  for  the  wavelets  and  packets  (equation 
21  and  27)  were  derived  from  the  ray  equations  of  geometrical  optics 
(equations  5  and  6) .  and  from  expressions  for  the  angular  frequencies 
of  the  wavelets  and  packets  in  a  moving  medium  (equations  13  and  14). 

In  geometrical  optics,  it  is  assumed  that  wave  properties  such  as 
wavelength  and  direction  change  little  over  a  distance  of  the  order 
of  a  wavelength. 

It  was  found  that  wave  properties  can  be  expressed  as  the  sum 
of  a  component  relative  to  the  current  and  a  component  due  to  advection 
by  the  current.  This  is  seen,  for  example,  in  equations  (13)  and  (14) 
where  the  angular  frequency  is  defined  with  respect  to  a  fixed  point. 
Further,  it  was  found  that  the  angular  frequency  relative  to  a  fixed 
point  is  constant  along  a  ray.  but  that  the  angular  frequency  measured 
relative  to  the  current  is  Doppler  shifted  by  the  current.  This 
latter  frequency  decreases  in  a  following  current  and  increases  in 
an  opposing  current.  Other  wave  properties  are  affected  in  a  similar 
manner.  In  particular,  the  ray  speed  P.  the  geometric  group  speed 
G' .  the  wavelet  speed  V*.  and  the  wavelength  increase  in  a  following 


current  and  decrease  In  an  opposing  current. 

Properties  of  wave  packet  rays  were  deduced  from  the  ray 
equations  (30)  and  (74).  From  these  equations,  it  was  found  that 
swell  waves  are  reflected  or  deflected  from  their  original  paths  by 
a  current  with  the  short  period  waves  being  affected  to  a  greater 
extent  than  the  long  period  waves.  Rays  were  found  to  bend  away 
from  (.toward)  the  current  normal  in  a  following  (opposing)  current. 

This  corresponds  to  saying  that  the  sign  of  the  curvature  of  a  ray 
is  the  same  as  the  sign  of  the  vortlcity  of  the  current.  The 
curvature  of  a  ray  was  also  found  to  increase  as  the  vorticity  of 
the  current  Increased  and  to  decrease  as  the  value  N  decreased. 
Further,  it  was  found  that  for  a  following  and  an  opposing  current 
the  change  in  the  direction  of  the  ray  is  greater  than  the  change 
in  the  direction  of  the  packet,  which  is  greater  than  the  change  in 
direction  of  the  wavelets.  If  u  «  G',  then  p  =  6.  Wavelets  and 
packets  were  found  to  remain  normal  to  the  current  contours  if  they 
entered  the  current  at  normal  incidence,  but  the  ray  was  found  to 
curve  due  to  advection  by  the  current.  It  was  found  that  In  a  parallel 
current  only  rays  in  a  following  current  are  totally  reflected,  and 
at  the  turning  point  the  ray,  packet,  and  current  directions  are 
parallel.  For  a  parallel,  opposing  current  rays  were  found  to  bend 
toward  the  current  normal  and  if  sufficient  refraction  occurred  rays 
would  become  perpendicular  to  the  current  contours.  For  curved 
currents,  it  was  found  that  total  reflection  could  occur  in  an 
opposing  current.  When  this  happens  the  ray  and  current  are  in 
opposite  directions. 


Snell's  law  equations  for  wavelets  (equation  37)  and  for 


packets  (equation  42)  were  deduced  from  the  fact  that  relative  to  a 
fixed  point  the  angular  frequencies  u  and  a  are  constant  along 
a  ray,  and  that  the  components  of  the  wave  numbers  k  and  m  parallel 
to  the  current  contours  are  also  constant.  When  the  transmitted  value 
of  the  angle  (6  or  y)  is  set  equal  to  90  degrees,  the  Snell's  law 
equation  can  be  rearranged  to  obtain  an  expression  for  the  critical 
angle  (wavelets,  equation  51  and  packets,  equation  53).  It  was  found 
that  critical  angles  are  dependent  upon  the  value  of  M  only.  Thus 
waves  of  various  periods  may  have  the  same  critical  angle  if  the 
current  values  are  sufficient  to  cause  the  value  of  M  to  be  the 
same.  For  example,  a  ray  with  an  Initial  angle  of  71  degrees  is 
reflected  if  M  is  0.046.  This  corresponds  to  a  14  second  period 
wave  and  a  current  speed  of  0.5  m/sec,  or  a  7  second  period  wave  and 
a  0.25  m/sec  current  speed.  This  relationship  is  shown  graphically 
in  Figures  2  and  3.  From  these  figures  it  was  found  that  the  critical 
angles  for  wavelets  are  greater  than  the  critical  angles  for  packets 
for  any  given  M  value.  Since  the  packets  and  rays  totally  reflect 
at  the  same  point,  whereas  the  wavelets  do  not,  the  values  of  8  and 
y  are  appreciably  different  following  a  turning  point.  This  would 
affect  the  geometric  group  speed  G'  causing  it  to  be  reduced  in  value. 
In  fact,  G'  will  vanish  if  6  and  y  differ  by  90  degrees.  Further, 
it  was  found  that  for  swell  waves  (7  to  20  seconds),  rays  with  Initial 
angles  less  than  44  degrees  are  not  reflected  by  open  ocean  currents 
since  the  current  speeds  are  not  sufficiently  large.  The  critical 
angles  computed  from  equations  (51)  and  (53)  are  for  currents  with 


parallel  contours.  It  was  found  that  the  critical  angle  for  waves 
which  enter  a  curved  current  from  the  inner  radius  side  are  larger 
than  the  critical  angles  determined  for  the  same  waves  in  a 
parallel  current,  and  that  the  critical  angle  for  waves  which  enter 
from  the  outer  radius  side  are  less  than  those  determined  for 
parallel  currents. 

Four  current  models  were  considered  in  this  study.  These 
models  were  patterned  either  after  the  Gulf  Stream  or  after  the 
Circumpolar  Current.  Two  of  the  models  had  parallel  speed  contours 
and  two  had  curved  speed  contours  (.Section  IV .A.).  Ray  trajectories 
computed  using  these  current  models  demonstrated  the  physical 
properties  determined  from  the  ray  equations.  For  the  parallel  current 
models.  It  was  found  that  rays  were  deflected  or  reflected  by  a 
following  current,  and  only  deflected  by  an  opposing  current.  Every 
ray  whldh  entered  an  opposing,  parallel  current  passed  through  the 
current,  but  was  deflected  some  distance  upstream  from  its  non- 
current  path.  Some  rays  in  a  following,  parallel  current  passed 
through  the  current,  and  others  were  turned  back  and  exited  the  current 
on  the  same  side  they  entered  the  current.  Similar  results  were 
obtained  from  the  curved  current  models  with  a  few  exceptions.  It 
was  found  that  some  rays  were  totally  reflected  by  an  opposing,  curved 
current.  Only  rays  which  entered  the  opposing,  curved  current  from 
the  outer  radius  side  could  possibly  be  totally  reflected.  It  was 
also  found  that  some  rays  were  focused  by  the  current  into  a  given 
area,  that  some  rays  exited  and  then  re-entered  the  current,  and  that 
some  rays  were  trapped  within  the  current. 
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These  results  for  the  parallel  and  curved  current  models  give 
some  Insight  as  to  what  could  be  expected  for  waves  in  the  actual 
Gulf  Stream  or  Circumpolar  Current.  All  waves  which  enter  the  Gulf 
Stream  from  the  northeast  would  be  deflected  some  distance  upstream 
from  their  non-current  path.  The  degree  of  deflection  would  depend 
upon  the  value  of  N  and  upon  the  curvature  of  the  current.  Waves 
which  enter  the  Gulf  Stream  from  the  southeast  may  pass  through  the 
current  or  may  be  totally  reflected  by  the  current.  Rays  with 
initial  angles  greater  than  about  44  degrees  may  he  totally  reflected. 
It  was  shown  that  the  critical  angle  increases  as  N  increases  and 
as  the  curvature  of  the  current  Increases.  Waves  which  enter  the 
Gulf  Stream  from  the  northwest  or  southwest  may  pass  through  the 
current  or  may  be  reflected  by  the  current.  The  greater  the  cur¬ 
vature  of  the  current  contours  the  more  likely  the  waves  which  enter 
from  this  side  (outer  radius  side)  of  the  current  will  be  reflected. 
These  results  show  that  the  Gulf  Stream  can  act  as  a  shelter  for  the 
coastline,  keeping  some  swell  waves  from  reaching  the  coast  from  the 
open  oceans .  However,  it  is  also  seen  that  some  waves  which  approach 
from  the  slope  water  side  will  be  turned  back  toward  the  coastline 
and  possibly  be  trapped  between  the  current  and  the  coastline. 

For  the  Circumpolar  Current,  it  would  be  expected  that  waves 
which  enter  from  the  Antarctic  side  would  pass  through  the  current 
with  some  deflection  from  their  non-current  paths.  Some  rays  may  be 
totally  reflected,  but  only  if  their  initial  angles  are  as  large  as 
about  75  degrees.  Waves  which  enter  the  Circumpolar  Current  from 
the  outer  radius  side  (from  the  north)  are  more  likely  to  be  reflected 
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by  Che  current.  These  rays  may  also  be  focused  by  Che  current  or 
exit  and  reenter  Che  current  on  the  Antarctic  side.  Thus,  the  Cir¬ 
cumpolar  current  may  significantly  alter  the  path  of  waves  which 
enter  it. 

Lastly,  the  wave  packet  refraction  method  developed  in  this 
study  was  compared  to  Arthur's  (1950)  monochromatic  refraction  method 
and  to  Kenyon's  (1971)  and  Teague's  (1974)  wave  packet  refraction 
method.  It  was  found  that  when  the  value  of  N  is  very  small 
(£  0.10)  and  when  the  initial  angles  of  the  rays  are  small  (<  30 
degrees),  that  the  various  wave  refraction  methods  give  nearly 
identical  trajectories.  However,  as  N  Increases  and  as  the 
initial  angle  increases  the  trajectories  become  considerably  different. 
The  differences  in  the  trajectories  are  due  primarily  to  the  differ¬ 
ences  in  0  and  y,  which  are  only  accounted  for  in  the  wave  packet 
refraction  method  presented  in  this  study. 

B.  Recommendations 

This  study  was  confined  to  deep  water  currents,  while  Arthur 
(1950)  considered  shallow  water  currents.  A  natural  extension  of 
this  study  would  be  to  investigate  wave-current  interactions  at 
intermediate  water  depths.  This  study  was  also  confined  to  con¬ 
sidering  ray  and  wave  directions  only.  Thus,  the  effect  of  the  cur¬ 
rent  on  wave  heights  was  not  considered.  Johnson  C1947)  and  Arthur 
(1950)  have  shown  that  wave  heights  are  altered  by  a  current.  Wave 
heights  may  be  chanced  h*  an  opposing  current  or  diminished  by  a 
following  current  Foi  sample,  Johnson  (1947)  discussed  waves  off 
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Che  coasC  of  Washington  which  were  observed  to  be  diminished  in 
height  by  a  current.  Arthur  (1950)  showed  that  wave  heights  could 
be  increased  such  that  breaking  occurred  within  the  current.  Thus, 
to  obtain  a  good  understanding  of  wave-current  interactions  wave 
heights  must  be  considered.  Another  area  which  needs  to  be 
investigated  is  the  effect  that  waves  have  on  the  currents.  Waves 
and  currents  probably  affect  each  other.  Teague  (1974)  mentioned 
the  possibility  that  waves  which  break  within  an  eddy  may  enhance  the 
decay  of  the  eddy.  Finally,  it  should  not  be  assumed  that  currents 
are  constant  to  the  depth  of  penetration  of  the  waves.  This  assump¬ 
tion  is  good  for  the  Circumpolar  Current,  since  it  is  a  deep  current, 
but  it  may  not  be  too  good  for  the  Gulf  Stream. 

It  is  highly  desirable  to  conduct  a  field  study  in  order  to 
test  the  refraction  theories  for  actual  ocean  currents.  The  Gulf 
Stream  would  be  a  good  current  for  investigating  wave-current  inter¬ 
actions,  since  its  current  speeds  are  large.  It  has  been  found  in 
this  study  that  large  current  speeds  show-up  the  differences  in  the 
various  refraction  methods  quite  well.  To  test  these  methods  good 
directional  wave  data  and  current  speed  profiles  would  be  required. 
Wave  data  should  be  taken  on  either  side  of  the  current  and  possibly 
at  a  couple  of  points  within  the  current.  Thus,  the  waves  which 
pass  through  the  current,  those  which  do  not,  and  the  amount  of 
deflection  from  their  original  paths  could  be  determined.  The  study 
would  be  enhanced  by  considering  waves  with  Initial  angles  in  the 
aid-range  (30  to  70  degrees),  since  rays  with  these  initial  angles 
differ  the  most  for  the  different  refraction  methods.  Wave  tank 
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experiments  could  also  be  conducted.  This  problem  lends  Itself 
quite  well  to  wave  tank  studies,  since  short  period  waves  are  affected 
the  most  by  a  current.  It  is  necessary  to  use  short  period  waves  in 
wave  tanks  in  order  to  have  sufficiently  short  wavelengths. 
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the  ray.  Three  equations  are  necessary,  since  after  refraction  the  wavelet  and 
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racket  directions  differ,  and  a  ray  is  not  generally  orthogonal  to  the  wave 
fronts  in  a  current.  These  equations  are  derived  and  applied  to  idealized 
current  models.  The  current  models  are  patterned  after  the  Gulf  Stream  and  the 
Circumpolar  Current.  Both  parallel  and  curved  current  contours  are  considered. 
From  this  study,  it  was  found  that  wave  packets  can  be  totally  reflected  or  de¬ 
flected  from  their  original  path  by  a  current  with  the  short  period  waves  being 
affected  to  a  greater  extent  than  the  long  period  waves.  The  curvature  of  a  ray 
was  found  to  have  the  same  sign  as  the  vorticity  of  the  current  and  to  increase 
as  the  vorticity  increased.  Rays  were  found  to  bend  away  from  the  current 
normal  in  a  following  current  and  toward  the  current  normal  in  an  opposing 
current  becoming  perpendicular  to  the  current  speed  contours  if  sufficient  re¬ 
fraction  occurred . \Wave  refraction  was  found  to  be  the  least  for  rays  with  small 
initial  angles  and  tor  currents  where  u  /G.  (maximum  current  speed  to  initial 
geometric  group  speep)  is  small.  Furthef*  it  was  found  that  the  change  in  the 
direction  of  a  ray  was  greater  than  the  change  in  the  direction  of  the  packet, 
which  was  greater  than  the  change  in  the  direction  of  the  wavelets.  The  ray 
angle  increased  in  a;  following  current ,  and  decreased  in  an  opposing  current . 
Critical  angles  were'  computed  for  parallel  current  speed  contours.  From  an 
analysis  of  these  angles,  it  was  found  that  wave  packets  totally  reflect  before 
the  wavelets,  and  at  the  same  point  as  do  the  rays.  Also,  critical  angles  were 
found  to  be  a  function  of  the  ratio  u  /G.  (current  speed  at  the  turning  point  to 
the  initial  geometric  group  speed),  aSd  to  increase  as  this  ratio  decreased. 
Further,  it  was  found  chat  rays  with  initial  angles  less  than  44  degrees  would 
not  be  totally  reflected  by  open  ocean  currents,  since  extremely  large  current 
speeds  would  be  required.  An  analysis  of  rays  in  currents  with  curved  contours 
indicated  that  some  rays  are  focused,  that  some  rays  exit  and  re-enter  the 
current,  and  that  some  rays  may  be  trapped  within  the  current.  These  are 
features  which  cannot  occur  with  parallel  currents.  Three  additional  refraction 
methods  were  considered  and  compared  to  the  method  presented  here.  They  were  a 
monochromatic  refraction  method  patterned  after  Arthur  (1950),  a  wave  packet 
refraction  method,  using  conventional  group  speed,  developed  by  Kenyon  (1971), 
and  a  vorticity  equation  used  by  Teague  (1974)  and  presented  by  Kenyon  (1971). 

It  was  found  that  there  was  very  good  agreement  in  the  results  of  all  the 
methods  for  small  initial  angles  and  small  current  speeds.  However,  differences 
between  the  results  obtained  by  the  different  methods  increased  greatly  for  mid¬ 
range  angles  (30-75  degreed  and  for  large  current  speeds. 
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